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ABSTRACT 


Mobiles are a communication facility connecting people from around the world. As the 
use of mobiles has increased and the technology advanced, analyzing mobile phone 
communication networks is very interesting. The data collected from the mobile 
companies and the technological development made in the area help understand human 
interactions pattern and the society. In the present paper, we used data records of calls 
and SMS from a mobile phone company in Albania to construct and analyze these two 
communication networks: i) the directed network (DN), by linking callers or SMS senders 
and call or SMS receivers with directed edges and, ii) the undirected mutual network 
(MN) considering only the pairs of users with reciprocal communication events. In 
both cases, we assigned a weight on each edge, based on the number of communications 
occurred for the corresponding pair of users during the observed time period. Based 
on a carefully chosen set of topological and weighted measures, we present a detailed 
analysis of the networks under consideration, by exploring the component size, degree, 
strength, and weight distributions, topological and weighted degree assortativity and 
clustering coefficient. The correlations between these quantities were inspected, 
providing a better insight into the local structure and interactions, and the overall 
topology of the networks. This is the first study that attempts to analyze and discover 
the topological and weighted properties of mobile phone communication networks 
with Albanian users. 

Keywords: mobile phone communication, network analysis, weighted networks 


1. INTRODUCTION 


Online social networks (Facebook, Twitter), E-mail, and mobile phones have 
transformed the way people communicate and interact with each other, providing 
connections with everyone all the time. The digital traces of social interactions, 
their frequencies and volume of communications generated by millions of people 
via all these communication channels offer a comprehensive picture of human 
habits. The management of such huge databases, network theory and statistical 
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tools have radically changed the scientific approach towards the understanding 
of our lives, the dynamic and performance of individuals, groups, organizations 
and societies (Watts, 2007; Lazer et al., 2009; Giles, 2012). 

An anonymized data set of calls and text messages records generated by the 
users of a cellular operator in Albania is here used to explore the Albanian 
mobile phone communication networks. We performed a detailed analysis of 
the networks under consideration, based on a broad set of measurements that 
researchers have proposed to express and characterize the topological and 
dynamical features of the real-world networks (Newman, 2003; Barrat et al., 
2004; Onnela et al., 2005; Boccaletti et al., 2006; Saramiaki et al., 2007). 


2. MATERIALS AND METHODS 


Data set The data consists of all call and text message (SMS) records among 
the subscribers of a cellular company which operates in Albania, during two 
nonconsecutive months or 61 days. This is a triple record ((i, 7), f), i.e., the user 
i (anonymized) calling or sending a SMS to the user j (anonymized) at date/ 
time f. There are 4,323,433 records, where 3,910,079 are voice calls and 413,354 
are SMS. 

Construction of networks We constructed the original directed network 
denoted as the DN. The nodes represent the active mobile phone users, and 
there is a directed edge from a node i to a node j if the user i has called or sent 
a SMS to the user j at least once during the observed period. We excluded the 
loops from the network. There are 102,872 nodes and 243,197 edges in DN. 

An undirected mutual network (MN) is obtained from the DN network as 
follows: each pair of edges (i, j) and (J, i) is substituted with an edge ij in the 
network MN, and all other edges are eliminated. The nodes with degree zero 
are removed. By this restriction, we only retain reciprocated ties and disregard 
the pairs of nodes connected by unidirectional edges, which mostly correspond 
to single communication events, suggesting that these pairs of users might not 
know personally each-other or the callers or SMS senders are more likely to 
be telemarketers or costumer service lines etc. Therefore, the mutual network 
MN should reflect a more realistic network of social interactions between the 
involved individuals. This filtering leads to the removal of about 44% of the 
total edges and 48% of the total nodes in the DN, remaining 53,453 nodes and 
68,055 edges in the MN. 

We assign a weight on each edge of both networks DN and MN to quantify 
the strength of a connection. In the DN, the weight w, of a directed edge (i, j) 
is the total number of calls and text messages initiated from the user i to the 
destination, the user j over the observed time period, whereas in the MN the 
weight w,, of the edge ij is the total number of communication events (calls or 
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SMS) occurred between the users i and j. 

Fitting distributions All the plots are in log-log scales unless otherwise 
specified. The scaling exponent of the power law distribution that best fit the 
component sizes 1s found using the Maximum Likelihood Estimation (MLE) as 
given in (Clauset et al., 2009). The best fitting parameters of the DPLN 
distribution are found with the MLE method, using the fitdistrplus package in R 
referring to the definitions and the estimation method proposed by Reed and 
Jorgensen (2004) and further information is provided in Noka (2014). 

The exponent of the power law If there is a power law dependence of 


some data y on x of the form: y = Cx“, the plot of the function y against x in 


doubly-logarithmic scales is a straight line with slope a. In the cases we treat 
here, we do as follows: We bucketize the x axes in logarithmic binning (Newman, 
2006). For each bin we calculate the mean of y values, and plot them versus 
the midpoints of the bins in log-log scales. Using a linear fit, we find the exponent 
of the power law, if there is such dependence. All of the exponents are fitted 
with R? > 0.95. 


3. RESULTS 


Connected components 

Both constructed networks are fragmented in a large number of connected 
components where the largest one contains the most of nodes as a giant 
component. In the directed network DN, 91,586 or 89% of the total nodes 
belong to its giant component, denoted as the g-DN. The giant component 
gcMN of the mutual network MN has 39,639 nodes or 74.2% of the total 
nodes. The mutual network MN is more fragmented than the DN, due to the 
removal of unidirectional edges. Table 1 summarizes the aformentioned 
information giving also the number of connected components and the number 
of nodes in the second largest component N,, . for both networks. 


Table 1. Sizes of the networks and their giant components 


The giant 
component Nstc 
| Network Nodes | Edges components Nodes | Edges | 


DN 102872 | 243197 4307 91586 | 233970] 21 


Connected 


MN 53453 | 68055 5002 39639 | 58718 26 


The component size distribution P(S) is the probability that a randomly chosen 
component in a network has S nodes. Figure 1A depicts the empirical component 
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size distributions plot for the networks DN and MN excluding their giant 
components. Both distributions exhibit power law distributions of the form 


P(S)~S~@, with exponents a,,,=3.52, a,,,=3.298. 


10" 5 
Component size S Distance | Distance | 


Fig. 1. Component size distributions (A). Number of nodes N (/) for several samples and 
their average versus the distance / in the gcDN (B) and gcMN (C). 


Snowball sampling is a method for obtaining samples of the large networks (Lee 
et. al., 2009). It starts from a randomly chosen node v, and includes all the nodes 
within a distance / (length of the shortest path) from the source node v. For a 
randomly selected set of source nodes, we obtain some samples using this method 
in each of the giant components gcDN and gcMN. By denoting as N (J) the number 
of nodes in a sample, we plot each of them (solid lines) and their average (dashed 
red line) as a function of the distance /=1,2,3,... in log-linear scales. The illustrated 
results are shown in the panels B and C of Figure 1 for the gcDN and gcMN, 
respectively. The dashed red lines show an exponential growth up to the distance 
[=6, followed by a slower increase rate to the total number of nodes of the respective 
networks (dashed black lines), indicating the presence of the small-world effect 
(Watts and Strogatz, 1998) in both giant components gcDN and gcMN. 


Fig. 2. A sample of the network gcMN, for / = 6. 
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A sample of the giant component gcMN, for /=6 is in figure 2 depicted. The 
color of nodes (yellow-green-blue-magenta-orange-lavender) corresponds to 
their distance from the source node colored in red. The lavender nodes dominate 
over all other nodes, whereas the neighbors of the red node are only 3 (yellow 
nodes). The width and the color (from yellow-green-cyan-blue to red) of the 
edges are drawn based on their weights. 

Basic properties of the nodes and edges. The basic characteristic of the 
nodes in a network is the degree distribution denoted as P(k), which gives the 
probability that a randomly chosen node has degree k. Panels A and B in figure 
3 show the degree distributions of the networks DN, MN and those of their 
giant components gcDN, gcMN, respectively. The degree distributions of 
directed networks DN, gcDN have a fatter tail than those of the mutual networks 
MN and gcMN. The maximum degree of the gcDN is k, = 1010, and the 
gcMN has k,,. = 77, which is a more “reasonable” number of contacts for a 
single individual. The average degree of the g¢-DN and gcMN are (k) = 5.1 
and (k) = 2.96 respectively, meaning that a user averagely communicates with 
about 5 persons, while he communicates reciprocally with about 3 of them. 

All these distributions can be well fitted by the Double Pareto Lognormal 


Distribution DPLN (@, f, V, T) with probability density function: 


_ a8 _ at? Inx— v— at? 
fix) a E texp (a a )o(==), 


ae inxs —v + Br? 
Pep (—pv +t )oe(oetFr) 


x >0,a,8,v,t > 0. 


@ is the cdf of a standard normal distribution and ®* = 1 — @ its ccdf. 

The estimated best fitting parameters are as follows: 

DN @ = 8.28, 6 = 0.1,v = 2.1,t = 0.76 and for its giant component gecDN 
@=83,8 =01,v=22,7=0.75. 

MN a@ = 9.4,8 = 0.1,v = 1.4, t = 0.65, and for the giant component g¢MN 
a= 9.17,6 =02,4=—15;1 = 03.45. 


Each network and its corresponding giant component present approximately 
identical degree distributions. From now on, we will focus on the analysis of 
the giant components gcDN and gcMN. 
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Fig. 3: Degree, strength, weight distributions fitted with DPLN distribution. 


The panel C of figure 3 shows the in-degree and out-degree distributions of 
the gcDN fitted by DPLN distributions with the corresponding parameters: 


a" =8665 6" = 029" = "= , 
aout = 7,9, B°¥F = 0.2,v°™ = 1.6, 7°“ = 0.79. It is observed that the nodes 


with very high degrees in the gcDN are mostly the contribution of their out- 
degree indicating that these hubs might be telemarketers or customer service 
numbers. 

The communication intensity of any pair of users reflects the strength of 
their relationship, and numerically is characterized by the weight of the 
corresponding edge in the networks under consideration. Not only the node 
degrees, but also the tie weights show a heavy-tailed distribution in both networks 
gcDN and gcMN, as seen in Figure 3D. This indicates a strong heterogeneity 
in the way people choose to distribute their communications across their social 
contacts. Most of edges have small weights (weak ties) and only a few of 
them have high weights (strong ties) (Granovetter, 1983; Onnela et. al., 2007a). 
Two connected users have reciprocally exchanged averagely 44 calls or SMS 
and maximally 5572 calls or SMS. The DPLN distribution fits well also the 
weight distributions in both networks gcDN and gcMN with the best fitting 


parameters as follows: for the gc-DN @ = 4.1,8 = 0.001,v = 2.7,7 = 1.42, 
and for the g¢MN @ = 4.1,6 = 0.3,v = 3.28,t = 1.47, 
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The strength of a node i, 5i = > jeN(i) Wij is an extended definition of the 


node degree. In the networks under consideration, the strength of a node is the 
total number of calls or text messages made or received by a certain user 
during the observed time period. This quantity measures the volume of 
communications that a user performs and his importance (Barrat et. al., 2004). 
Strength distributions for the networks gcDN and gcMN are plotted in Figure 
3E. For the directed network gcDN, we distinguish the incoming and outgoing 
strength of nodes by summing up the weights of incoming and outgoing edges 
respectively, and their distributions are plotted in Figure 3F. Obviously, there is 
a large variability of the volume of communications that people perform through 
their mobile phones in each direction too. Although the node strength is 
averagely 81.3 (1.3 events per day) in the gcDN and 130.6 (2.1 events per 
day) in the gcMN, the maximum value is 12402 (203.3 events per day) and 
6651 (109 events per day) respectively. All the strength distributions can be 
fitted very well with DPLN distribution with the following parameters: for the 


gcMN a= 7.56,B = 0.3,v = 5.3,t = 1.05, for the gcDN, 
a= 7.2,8 = 0.05,v = 5.5,t = 0.9905, for the incoming strength in the gcDN, 
a= 7.32,68 = 0.01,v =4.97,1=0.9, and for the outgoing strength 
a= 7.1,68 = 0.35,v = 4.1,7 = 1.45. 

Similar pictures of the degree, strength and weight distributions observed in 
both networks gcDN and gcMN suggest the existence of mutual dependences 
between these quantities, leading us to a more detailed analysis below. 

Degree, strength and edge weight correlations Degree correlation is an 
important parameter to better understand the topology of the networks. We 
want to investigate if there is a relation between the number of incoming and 
outgoing neighbors for a single node in the directed network gcDN. The scatter 
plots of in-degrees k-’" and out-degrees k-*" per node are in Figure 4A depicted. 
The blue dots are the medians by logarithmic binning (see Materials and 
Methods). It is observed the power law dependence k°“* ~(k'")® with 
a=1.045. This means that the more popular individuals (those receiving many 
calls or text messages) tend to contact more people. 

A measure of degree-degree correlation is the average nearest neighbors 


degree of a node i: 
ae 
nn,i k; jen(i) j 


where N(i) is the neighborhood of node i. By calculating the average of the 


Knni values over the nodes with degree k, one can get the average degree of 
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Fig. 4. In-Out degree correlation in the gcDN (A). Degree and weighted degree correlation 
for the gcMN (B) and the gcMN (C). 


the nearest neighbors denoted as (Kn»| &}. This quantity can be also expressed 


as (Knn| kK) = YLerk’P(k'|k), where P(k’|k) is the conditional probability that 
a node with degree k’ is connected to a node with degree k. In absence of 
degree correlations, P(k’|k) does not depend on k, and neither does , hence 
(Pastor-Satorras et al. 2001). On the contrary, if increases as a function of k, 
high-degree nodes have a lager probability to be connected to other high-degree 
nodes. In this case, the network is said to show degree assortative mixing 
(Newman, 2002). Barrat et al. (2004) has introduced the weighted average 
nearest neighbors degree of a node i in the weighted networks: . By averaging 
the values over the nodes with degree k, the behavior of the function marks 
the weighted assortative or disassortative properties, considering the interactions 
among system’s elements. In Figure 4B and C, we plot and as functions of 
degree k for the networks gcMN and gcDN respectively. For the gecDN, we 
do not consider the direction of edges. In the mutual network gcMN, both 
functions and have an increasing trend with k up to the k = 27 and k = 34, 
respectively. The nodes with degree larger than k = 27 constitute only 0.08% 
of the total nodes. Therefore, the mutual network gcMN shows the degree 
(weighted) assortative properties, indicating that the users preferentially establish 
their relationships with individuals similar to themselves. This human tendency, 
known also as affinity or homophily emerges also with respect to a variety of 
attributes such as health or habits (Christakis and Fowler, 2007), tastes and 
interests (Lewis et al., 2008), age or race (Mollica et al., 2003). In the directed 
network gcDN, the function shows an increasing trend along a wider range of 
the k values than the , but the nodes with degree within this range constitute 
only 18.3% of the total nodes of the network. Thus, gcDN exhibits mostly 
disassortative mixing pattern. 

We can also study the correlations between nodes strength. The average 
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nearest neighbors strength of a node / is defined as . By averaging values over 
the nodes with strength s, the average strength of nearest neighbors is obtained. 
We plot as a function of the strength s in Figure 5A for the networks gcDN 
and gcMN. Both curves present qualitatively similar pictures. The function 
shows an independence for small s values, followed by a slightly increase with 
the growth of s and a more obvious growth from s~800 onwards. The nodes 
with strength greater than 800 establish only 2.5% of the total nodes in the 
gcMN, and 1.1% in the gcDN. Furthermore, the strength of such nodes is 
dominated by the weight of a single edge jj, thus, s.~w,~s,,. These results 
suggest that in both networks, the strengths of two adjacent nodes are mostly 
uncorrelated. 

The similarity of the strength distribution P(s) with the degree distribution 
P(k) in both networks, as seen above, leads us to analyze the relations between 
the node strength and the node degree in the networks gcDN and gcMN. If 
the weights of the edges incident to a node would not correlate with its degree, 
then with a = 1, because , where is the average weight of the network. By 
plotting the average node strength as a function of degree shown in Figure 5B, 
we find that the node strength follow a power law growth up to k = 45 in the 
gcMN, with , and for larger k values the strength saturates. In the gcDN, the 
power law dependence exists with , up to k = 25, and for larger k values we 
observe a clear deviation from the linear approximation. These observations 
indicate the existence of strength-degree correlations and suggest that the 
individuals with many contacts averagely dedicate less calls or SMS per contact 
than those with few connections. The node strength conditional on in/out degree 
for the gcDN presents similar results (Figure 5C; the green dots are shifted 
vertically by a factor of 10 for better visibility). For small values, both curves 
exhibit power laws, where ai = 1.15 and a = 0.94, and deviations are 
observed for greater values of k. 

In addition to the strength-degree correlations, we plot the average strength 
product conditional on degree product , shown in Figure 5D for both networks 
gcDN, gcMN. In the case of absence of correlations, we would expect . Then 
with B = 1. On the contrary, it is found that and . 

Similarly, we analyze the average weight conditional on degree and strength 
product, and respectively. The corresponding results are presented in Figure 
5E and F. In the mutual network gcMN, the exponent of the power law of the 
form , is , which means that the strength of a tie (edge weight) between two 
individuals does not depend on the absolute number of friends that each of 
them has. Despite the slight increase of observed for small values , practically 
the network gcDN also shows an independence of the edge weight versus the 
degree product. In contrast, the rapidly increasing curves in Figure 5F indicate 
the presence of correlations between the edge weight and the strength product 
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Fig. 5: Average nearest neighbor strength (A). Average strength conditional on degree (B). 
Average out/in strength conditional on degree in the gcDN (C). Strength product conditional on 
degree product (D). Edge weight conditional on degree product (E) and on strength product (F). 


of its end nodes in both networks gcDN, gcMN with roughly power law 
dependences of the form . The scaling exponents are . 

Clustering coefficient. A local property of the network structure is the 
tendency of the neighbors of a node i to be interconnected, which is characterized 
by the (topological) clustering coefficient , where is the number of triangles 
around the node i. (Watts and Strogatz, 1998). The average clustering coefficient 
as a function of degree is presented in Figure 6A for the networks gcMN and 
gcDN. For the gcDN we do not consider the edge directions. Both curves 
show power law decay of the form for the degree values up to =40, which 
cover almost all the nodes. The power law dependence of the clustering 
coefficient around a node on its degree is observed in many empirical networks 
with different scaling exponents (Newman, 2003). 

In Figure 6B, we present the average weighted clustering coefficient as a 
function of the strength (Saramaki, 2007) for both networks gcDN, gcMN. We 
observe an increase of the triangles’ intensity around the nodes with the increase 
of their strength. But, it seems that the behavior of the high-strength nodes, as 
seen above for the , is different from that of the others. 
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Fig.6. Topological clustering coefficient as a function of degree (A). Weighted clustering 
coefficient as a function of strength (B). 


4. CONCLUSIONS 


We have constructed two networks from mobile phone calls and SMS records 
of an operator in Albania. The weights of edges used to measure the strength 
of ties as the total number of communications during the observed time period 
have offered a good opportunity in exploring the mobile phone users’ interactions 
and their relations with the topological properties of the networks. 

Both networks present a large heterogeneity in the node degree and strength, 
and in the weight of ties which is implied by their distributions. The hubs (nodes 
with very high degree) in the directed network are more than those in the 
mutual network. This suggests that the reciprocal communication criterion for 
the construction of the mutual network eliminates such abnormal nodes, most 
likely associated with costumer service numbers or telemarketers, providing a 
better representation of the real social interactions. 

We have found that people have the tendency to form groups with others 
that are similar to themselves by the number of contacts, while most of them 
disregard the volume of communications performed by their connections. By 
investigating the relationship between the communications that people perform 
and the number of their social contacts, we have found that people with large 
social circles need to perform more communications. However, the very large 
number of connections does not mean a proportional increase of the volume of 
communications. The individuals with many contacts dedicate averagely less 
communications to each of them than those with few connections. In addition, 
there exists a strong diversity in the way people schedule and manage the 
attention to their social relationships. The strength of a tie between two individuals 
does not depend on the absolute number of their social contacts, but it is positively 
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correlated with the number of communications they perform. 

Albanian mobile phone communication networks share many common features 
with other similar networks studied before. Specifically, our results for the 
mutual network are in line with those reported by Onnela et. al, (2007b), who 
analyzed a large mutual network based on the data from an operator of an 
undisclosed European country with all mobile phone calls between 7 million 
people during a period of 18 weeks. 

Li, et. al.( 2014) studied both directed and mutual networks, similar with 
those we analyzed in this paper, using a dataset with detailed call records of 
more than nine million different mobile phone numbers from one mobile operator 
in Shanghai. The observations for directed networks are qualitatively similar. 
They compared their results of the mutual network with the aforesaid network 
of the European dataset and observed that both networks share many common 
topological features, but also several discrepancies. The mutual network of 
Shanghai dataset has a different broadness of degree distribution and it is not 
observed a dependence of the topological clustering coefficient on the node 
degree as observed in the European and Albanian network. 

Our presented study is the first attempt in exploring and analyzing the topology 
and interactions’ features of mobile phone communication networks with 
Albanian users. 
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ABSTRACT 


Bifurcation analysis is a key tool for the analysis of dynamic systems in general and 
nonlinear systems in particular. If a parameter of a dynamic system changes, qualitative 
behavior of the system will change. New fixed points might emerge or current fixed 
points might disappear, or their stability properties may change. These sorts of 
qualitative changes in a dynamic system are called bifurcations, and the values of the 
parameters at which these changes take place are called bifurcation points. Bifurcation 
analysis also enables us to qualitatively estimate the behavior of trajectories without 
utilizing the analytical solution of the underlying differential equations. Most common 
types of bifurcations encountered in dynamical systems are saddle-node, transcritical 
and pitchfork bifurcation. In the present paper some of the types of bifurcations in 
one-dimensional dynamical systems and application of these bifurcations in population 
dynamics are discussed. The fish available for human consumption comes either from 
the ocean or the sea. On one hand, the natural supply cannot satisfy the human needs. 
On the other hand, cost of fish harvesting is increasing, making aquaculture an important 
source for fish supply. As Albania is rich in water resources, the market demand for 
fish could be easily met. Aquaculture production is playing an increasing role in 
satisfying the demand for human consumption of fish and fishery products. Once 
aquaculture production is increased, the endangered species would be safe, because 
interventions in the rearing process to enhance production are involved. The logistic 
growth model is used for population growth of fish, and the two following harvesting 
strategies are considered: i) constant and, ii) proportional harvesting. In each process, 
the optimal amount of fish harvested to protect the population from extinction was 
estimated. The results reported that harvesting in amount or rate higher than the 
bifurcation point leads to the extinction of population. The results are a means to 
address the growth of fish population and reduce of repopulation’s costs. 
Keywords: fixed point, stability, saddle- node, harvesting, logistic growth model 
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1. INTRODUCTION 


Bifurcation analysis of dynamical systems helps understand the dynamic 
behavior of the systems and find out possible dependence of the system’s 
behavior upon parameter values by reporting on the occurrence and changes 
in stability of fixed points, and helping to model these changes and transitions 
from stable to unstable case or vice versa as some parameters change. In 
addition, bifurcation analysis enables us to qualitatively estimate the behavior 
of trajectories without utilizing the solution of the underlying differential or 
difference equations. This is achieved by numerically approximating equilibrium 
solutions and their stability, even for problems that do not have analytic solutions. 

Bifurcation is the change in number of fixed points or periodic orbits, or in 
their stability properties, as a parameter varies. A bifurcation point is the value 
of parameter at which the change occurs. Bifurcations are classified according 
to how stability is lost. All types of bifurcations here reported are local 
bifurcations, because only the behavior of a dynamical system in the 
neighborhood of a single fixed point is affected. 

Population dynamics is of great interest for the commercial harvesting industry 
and biology, ecology and economics etc. Mathematical models have been used 
widely to estimate the population dynamics of animals and humans. In recent 
years, the use of mathematical models has been extended to agriculture sector 
especially in cattle farming to ensure continuous and optimum supply. The logistic 
growth model in term of harvesting has been used to study the fishery farming 
(Laham et al., 2012). The most important for successful management of 
harvested populations is that harvesting strategies are sustainable, not leading 
to instabilities or extinctions and produces great results for the year with little 
variation between the years (Aanes et al., 2002). Therefore, it can satisfy the 
market demand throughout the year. 

The fish available for human consumption comes either from the ocean or 
the sea. Unfortunately, the natural supply cannot satisfy the human needs. On 
the other hand, cost of fish harvesting is increasing, making aquaculture an 
important source for fish supply. Aquaculture production is playing an increasing 
role in satisfying demand for human consumption of fish and fishery products. 
Once aquaculture production is increased, the endangered species would be 
safe, as interventions in the rearing process to enhance production are involved. 

As Albania is rich in water resources, the market demand for fish could be 
easily met. In recent years, fish farming has developed in Albania and aquaculture 
production constituted 20.6% of total fish caught in 2013, around 14% in 2010 
and around 1% in 2001 (http://www. instat.gov.al/media/154623/tb19.xIsx). 

The use of mathematical models in fishery harvesting helps the aquaculturists 
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to estimate the fish population for a given period. In addition, the models provide 
effective solutions with regard to the consumer’s demand. 

In the present paper, the saddle-node, transcritical and pitchfork bifurcations 
that occur in dynamical systems is discussed and the bifurcation analysis to harvest 
models is applied using the logistic growth model and the two following harvesting 
strategies: 1) constant harvesting and, ii) proportional harvesting; to estimate the 
optimum quantity of fish for harvesting that can ensure the continuous supply for 
each strategy and, to compare the results obtained between strategies. 


2. ONE DIMENSIONAL BIFURCATION 


Consider a dynamical system that depends on one parameter, h: 


X= f(x,h) x(0) = Xo (1) 


When a change in a parameter results in a qualitative change in the dynamics 
of a nonlinear system, the system is said to have gone through a bifurcation. 
Bifurcations are classified by the way in which the fixed points of the function 
jf change their number, location, form, and stability. 

This paper discusses three common types of bifurcations: saddle-node, 
transcritical and pitchfork bifurcations (Strogatz, 1994; Lorenz, 1993; Shone, 
2002; Teschl, 2004). 

In a saddle-node bifurcation, as the bifurcation parameter passes through the 
bifurcation point, two fixed points disappear, so there are no fixed points 
afterwards. Before they disappear, one of the two fixed points is stable and the 
other is unstable. 

If the differential equation (1) satisfies the equations: 


F(X») =0 
of _ 
ax (X54) =0 


2 


Of 
Bye Cor Mo) #0 


O 
aa Fash) #0 


then this equation has a saddle node bifurcation at h = h,. 
The key example of a saddle node bifurcation is provided by the differential 
equation: 
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x=h-x° (2) 


Clearly, we have f (0,0) =0 and = (0,0) = 0. Fixed points of this equation 
x 
are given by h= x’, corresponding to a parabola in the xh plane. For h < 0, the 
equation has no fixed points and x is always negative. For h > 0, there exist 
two fixed points, one stable and one unstable. 

In dynamical systems, a bifurcation diagram shows the possible long-term 
values (fixed points or periodic orbits) of a system as a function of a bifurcation 
parameter in the system. The stable fixed points are represented by a green 
line and unstable ones by a red line. The bifurcation diagram for the saddle- 
node bifurcation indicates again that the two fixed points collide and annihilate 
at h = 0 and there are no fixed point for h < 0 (Figure 1). 


Fig.1:Bifurcation diagram for the saddle-node bifurcation. 


In a transcritical bifurcation, there are two fixed points, one stable and one 
unstable. When the bifurcation point is passed, there is an exchange of stability; 
the unstable fixed point becomes stable and the stable one becomes unstable. 

If the differential equation (1) satisfies the equations: 


f (XN) =9 


0 
© ty) =0 
Ox 
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Of 
ap 7070) =0 


af 
Bye 607M) x 0, 


Of 
and axah (XM) #0, 
then this equation has a transcritical bifurcation at h = h,. 
The simplest example of a transcritical bifurcation is provided by the following 
differential equation: 


x=hx-x° (3) 


Clearly, we have f(0,0) =O and J, 0) =0. The fixed points of this 
x 
equation are x = 0 and x = A, corresponding to two lines in the xh plane. For h 
< 0, the equation has two fixed points, x = 0 is stable and x = his unstable. The 
two fixed points collapse at h = 0 and for h > 0, x = 0 is unstable and x = h is 
stable. An exchange of stability has occurred at h = 0. 
The bifurcation diagram for the transcritical bifurcation is depicted in figure 2. 
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Fig. 2: Bifurcation diagram for the transcritical bifurcation. 


In a pitchfork bifurcation, there are two stable fixed points separated by an 
unstable one. Once bifurcation point is passed, only one stable fixed point 
remains. This type of bifurcation is called supercritical pitchfork bifurcation. 
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There is also a subcritical pitchfork bifurcation, in which there are two unstable 
fixed points separated by a stable one, until the bifurcation point is passed. 
Then there is only one unstable fixed point. 

If the differential equation (1) satisfies the conditions: 


Ff (%5'y) =9 
F(,.h,)=0 
Fs,.I)=0 
SL Oa.) #0 
SE oh) #0 


3 
and ashy) #0, 


then this equation has a pitchfork bifurcation at h = h,. 
As an example of a pitchfork bifurcation consider the following differential 
equation: 


x=hx-x° (4) 


Clearly, we have f(0,0)=0 and J, 0) =0. The fixed points of this 
x 
equation are given by x =0 and h =x’, corresponding to a line and a parabola 
in the xh plane. For negative values of h, the equation has one stable fixed 
point at x = 0. For h > 0, x = 0 is a unstable fixed point and the two new fixed 
points have been created at h = 0, one in each branch of the parabola h = x’. 
The bifurcation diagram for the supercritic pitchfork bifurcation is depicted 
in figure 3. 
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Fig. 3: Bifurcation diagram for the supercritic pitchfork bifurcation. 


Next section discusses the application of these types of bifurcation to fishery 
management. 


3. APPLICATION OF BIFURCATIONS: 
FISHERY MATHEMATICAL MODELS 


To illustrate the bifurcation analysis in a real world situation we are going to 
analyze some simple fishery management models. Consider the logistic growth 
equation to model a fish population in the absence of fishing: 


dP P 
ae 6 
ated 9 


where P is the size of the population, ris the growth rate due to reproduction 
and M is the carrying capacity of the environment. Now considering modeling 
the population and harvesting some of the population using some common 
harvesting strategies: constant harvesting and proportional harvesting. Constant 
harvesting is where a fixed numbers of fish were removed each year, while 
periodic harvesting is usually thought of a sequence of periodic closure and 
openings of different fishing grounds (Idels and Wang, 2008; Laham ef al., 
2012). In proportional harvesting, the quantity harvested is proportional to the 
population. Harvesting has been considered a factor of stabilization, 
destabilization, and improvement of mean population levels, induced fluctuations, 
and control of non-native predators (Michel, 2007). 


94 Valentina SHEHU, Alma SPAHO 2015/Vol. XX (2) 


Constant harvesting 


One of the simplest methods is to establish a set limit for harvesting. We 
assume that the dynamics of the population satisfies the logistic growth model 
(5) and that a constant harvesting, h, is added for removing a constant number 
of the fish over a given time interval. 

The mathematical model becomes: 


dP P 
~~ =rP|1-— |-h = 
res =| PO) =p, 6) 


The solution of this ordinary differential equation obtained using MAPLE is 
difficult to analyse, so we turn to a geometric analysis of the model. 
The fixed points, P*, are simply the solution of the equation rP*(1- P’/M) = h. 


« | 4hM 
The model (6) has two roots: P,, = sla +,|/M*— a) . The number of 
r 


fixed points depends on the sign of the term M? - 4hM/r. There are: two fixed 
points if M? - 4hM/r > 0 or equivalently 0 < h < rM/4; one fixed point if M? - 
4hM/r = 0 or equivalently h = rM/4; and no fixed point if M? - 44M/r < 0 or 
equivalently h > rM/4. 

To show the above results, conditioned by the unavailability of any real and 
reliable data, consider a lake or pond with a rare fish and suppose that the rate 
of population growth is r=0.7, carrying capacity is M = 1000. The fixed points 
for model (5), where h = 0, are P* = 0 and P* = M = 1000. 


The bifurcation diagram for the constant harvesting logistic model is depicted 
in figure 4. 
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Fig. 4: Bifurcation diagram for the constant harvesting. 


JNTS Bifurcation analysis as a useful tool to fishery management 25 


As the harvesting increases, the two fixed point move closer to each other 
with the lower fixed point remaining unstable and the upper fixed point remaining 
stable. As h moves toward the maximum growth rate of the logistic growth 
equation, the two fixed points coalesce at P* = M/2 = 500. When h > rMA = 
175, there are no fixed point, and the model shows that the population always 
goes extinct. This model shows a classic example of a saddle node bifurcation. 
The bifurcation values is at h = rM/4 = 175. 

The purpose of simple mathematical models applied to complex problems is 
to offer some insight. Here, the results reported that overfishing (in the model 
h > rM/4) during one year can potentially result in a sudden collapse of the fish 
catch in subsequent years, so that governments need to be particularly cautious 
when contemplating increases in fishing quotas. 


Proportional harvesting 

Another common form of harvesting is when one puts in a constant effort to 
harvest. In this case, the quantity harvested is proportional to the population. 
Thus, the mathematical model is written: 


dP P 
—=rP|1—-— |-hP = 


where again r is the growth rate, M is the carrying capacity with no harvesting 
and now / is the proportional rate of harvesting. 

Algebraic solution is complex and harder to interpret, thus we again turn to 
the geometric analysis of the model. The fixed points of (7) are the solution of 


* 


oo P ok * ra 
the equation: rP ( 4 =hP ,thatis,and P = MOSH, 


r 


. The extinction 


fixed point, P* = 0, is unstable for values of h < r. As h increases, the larger 
equilibrium (carrying capacity, M) shrink, but it remains stable for h < r. Forh 
= 0, P*=M= 1000. 

The bifurcation diagram for proportional harvesting logistic model is depicted 
in figure 5. 
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A 


Fig. 5: Bifurcation diagram for proportional harvesting. 


As the harvesting increases, the nontrivial fixed point move closer to extinction 
fixed point. As h moves toward the growth rate (0.7), the nonzero fixed point 
fades to zero, which implies there is extinction. This is comprehensible due to 
the harvesting rate approaching the growth rate. When h > r (= 0.7), the rate 
of harvesting exceeds the reproduction rate and extinction necessarily follows. 
This model shows a classic example of a transcritical bifurcation. The 
bifurcation point is h = r (= 0.7). 

The results reported that overfishing (h > r) in one year can potentially extinct 
the fish in the lake or pond. Consequently, the bodies involved in the area need 
to be particularly cautious when contemplating fishing quotas. 


4. CONCLUSIONS 


Differential equations containing parameters with approximately known values 
are commonly involved in practical applications. Consequently, investigating 
the solutions’ behavior and their dependence on the parameters is very important, 
as a slight variation in a parameter might have significant impact on the solution. 

The three common types of bifurcations here reported are: i) saddle-node 
bifurcation, it) supercritical pitchfork bifurcation and, 1ii) transcritical bifurcation. 
In each case, the bifurcation value is h = 0. These types of bifurcations also 
occur in higher-dimensional dynamical systems. 

Sustainable harvesting strategies and successful management of harvested 
populations help the fishery industry to grow both qualitatively and quantitatively. 
However, some little variation in fish population between the years could still 


JNTS Bifurcation analysis as a useful tool to fishery management 27 


be reported. Regarding the constant harvesting strategy, the fish population 
does not have enough time to recover if the constant harvesting is greater than 
the bifurcation value h = rM/4. Regarding the proportional harvesting strategy, 
the fish population will extinct if the proportional rate of harvesting is greater 
than the growth rate of the population, h > r. 

The results are a means to address the growth of fish population and reduce 
of fish repopulation’s costs. In addition, the annual overfishing has unavoidably 
impacts fishing in the subsequent years. Here, an appropriate national policy 
on fishing quotas would be crucial. 

Appropriate fishery harvesting strategies would satisfy the market demand. 
Commercializing the aquaculture could be a good alternative. 

However, some limitations could be here encountered: 1) theoretical results, 
ii) unreal parameters of logistic model for fish harvesting is reported. 

The application of bifurcation analysis to fish harvesting in real data is a 
means to address further researches in the area. In addition, other practical 
applications of bifurcations such as one-dimensional or higher dimensional 
dynamical systems must further. 
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ABSTRACT 


Cluster ions and aerosols are atmospheric particles, which participate and have a great 
impact on several global processes. On other hand, both these particles interact with 
each other via recombination and attachment processes. Monitoring the variation of 
their concentrations gives valuable information about their interactions and the impact 
on atmospheric processes that these particles take part. The present paper provides 
information about the variation of particle number concentrations of cluster ions and 
aerosol particles of sub-micrometric and super-micrometric size. Investigation was 
carried out from 2009-2011 at a site located in the Adriatic seashore and the results 
reported high values of cluster ions and aerosols. The source is anthropogenic 
activities, sea salt and other particles of long-range transport. 

Keywords, Jon concentrations, aerosol concentrations, seashore site, anthropogenic 
activities 


1. INTRODUCTION 


The major part of the lower troposphere consists of atmospheric particles 
such as cluster ions and aerosol. The formation and growth of atmospheric 
aerosols depend on several steps, namely nucleation, initial steps of growth 
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and subsequent — mainly condensational — growth. The main source is 
anthropogenic activity, but they are also formed naturally in seas, deserts, and 
forests. Recent data report that production of new aerosol particles by nucleation 
and subsequent growth by condensation of organic vapours emitted by trees is 
a frequent phenomenon that occurs in most atmospheric environments (Kulmala 
et al., 2004). 

Cluster ions differ from aerosol particles primarily by their smaller size or 
greater electrical mobility (Horrak, 2001). Other differences between these 
two groups relate to their electric charge, spatial distribution, vertical profiles, 
size distribution functions, etc. Atmospheric ions, according to their electrical 
mobilities may be divided into three main categories (table 1) (Horrak et al., 
2003). 


Table 1. Categories of atmospheric ions 


Diameter Range Mobility Range, 
nm em’ V's 7 
Small ions 0.36 - 1.6 1.3-0.5 
Intermediate ions 1.6-7.4 0.5 - 0.034 
Large ions 74 - 80 0.034 - 0.00041 


Aerosol particles, according to their size, could be divided into the following 
groups: 1) ultrafine (3 — 90 nm), ii) fine (90 — 1000 nm) and, iii) coarse (> 1000 
nm). Ultrafine mode may be divided also into two other modes; nucleation and 
Atiken (on threshold 20 nm). Unlike atmospheric ions, aerosol particles can be 
neutral or single/multi charged (Mandija, 2011-a). Table 1 reports that two 
large categories of atmospheric ions (intermediate and large ions) are charged 
ultrafine aerosol particles. Thus, the traditional border between cluster ions 
and charged aerosol particles is 1.6 nm (Tammet, 1995). 

Both atmospheric ions and aerosols play important roles in atmospheric 
processes (Flagan, 1998). Air ions make the air fresh, helping persons who 
suffer from asthma or other respiratory problems. Meanwhile atmospheric 
aerosols, “swallowing” air ions, play an inverse role on environmental concern 
(Rose et al., 2013). 

Concentrations of atmospheric ions and aerosols are governed by balance 
equation. This equation in its simplest form is presented by eq.1 (Salm and 
Tamm, 2011). 

dn 


7 4 ts — BunN, — Bn. — Bola. +2,) (1) 


where: n.,, is the concentrations of ions of different polarities 
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N,, ) is the concentrations of charged and neutral aerosols 


a and f are their respective recombination and attachment coefficients. 


Larger ions can swallow more cluster ions, but because of their minor 
concentrations, there are small aerosols that influence more of the ion 
concentrations. 

The present investigation was carried out at Velipoja beach which is a tourist 
location on the Adriatic seashore. Air studies are of great importance, as they 
help evaluate air quality (Mandija et al., 2010; Mandija, 2011-b). 


2. MATERIAL AND METHODS 


Measurement of atmospheric ion and aerosol concentrations was carried out 
at Velipoja, a tourist location on the Adriatic seashore in Albania. Figure 1 
depicts site location 200 m from the seashore, in the middle of a residential 
area. 

Investigation was carried out in spring, summer and autumn (April-October). 
In summer, anthropogenic activity intensifies along with the population, as several 
tourist-commercial activities take place. In spring and autumn, anthropogenic 
activities are very rare. 

Measurements were carried from 2009 to 2011. As only measurement results 
regarding fair-weather conditions were analyzed, the impact of precipitations 
and strong winds on atmospheric particle concentrations were excluded. 


Fig.1 Monitoring site, Velipoja Beach, Albania (from Google Maps). 


The Air Ion Counters IC-1000 and Andes ITC-201A were used to measure 
the concentrations of atmospheric cluster ions (d < 1.6 nm). Operation principle 
of this instrument is the appliance of electric field between three parallel 
conductive plates. The middle plate collects air ions which enter into the intern 
area of the AIC. Atmospheric aerosol concentration was measured using a 
Handheld Particle Counter, model Handilaz mini 301. This instrument measures 
aerosol concentrations divided into three channels; in the size range over 0.3 
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um. Thresholds between channels varied from 0.5 and 5 um. Aerosols on the 
first channel belong to fine mode, whilst those on the third channel belong to 
coarse group. Aerosols on second channel belonged both modes, but more to 
the fine mode (because their size distribution). Operation principle of this 
instrument is laser light scattering. 


3. RESULTS 


3.1 Ion concentrations 

Table 2 reports the average values of meteorological parameters during the 
monitoring process and temperature, relative humidity, atmospheric pressure 
and wind speed had normal values throughout the monitoring process. Usually, 
in summer temperature is higher and wind speed has lower values 


Table 2. Average values of meteorological parameters 


Summer Spring and Autumn 
Temperature (°C) 28.6 19.6 
Relative humidity (%) 62.3 57.3 
Atmospheric pressure (hPa) 1205.1 1027.3 
Wind speed (km/h) 1.2 3.7 


Table 3 reports the average values and the standard deviations of ion 
concentrations of both polarities. 


Table 3. Seasonal-averaged ion concentrations 


Summer Spring and Autumn 
Positive (cm™) 1752.0+843.0 1145.8+119.6 
Negative (cm”) 1550.04665.1 487.5+100.3 
Total (cm) 3302.0+1482.8 1633.3+182.6 
Unipolarity factor 1.1+0.2 2.5+0.7 


The average values of the concentrations of both positive and negative ions 
in summer were higher (positive 1.5 and negative 3.1 times) than in spring and 
autumn. Since the anthropogenic activities, like traffic, heating-cooking, etc. 
are extremely higher in summer due to the tourism air ion concentrations are 
higher (Yu, 2002). In addition to the anthropogenic activity, vegetation in the 
surrounding forest has a great impact (Jayaratne et al., 2011). This fact implies 
the hypothesis that the principal sources of air ions in Velipoja beach are 
anthropogenic activities and vegetation. Moreover, the dry summer season also 
enhances radon gas emanation from the ground surface, increasing its 
concentration in the regions near the ground (Hosoda et al., 2007). Both these 
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factors have a renewed presence in this region during summer seasons. 

Table 3 and 4 report that the meteorological parameters had a great impact 
on the ion concentrations. As the temperature is higher in summer, ion 
concentrations are higher. So, we can suppose that temperature influences the 
increase of ion concentration (Luts, 1998) indirectly due to the augmentation 
of radon exhalation (Porstendorfer, 1994; Hopke and Ramamurthi, 1998; Katoka 
etal., 1998). Regarding the wind, its speed affects directly the ion concentration 
via dispersion. These hypotheses are confirmed by the correlation coefficients 
between ion concentrations, temperature and wind speed. Table 4 presents the 
Pearson correlation coefficients among the ion concentrations and the principal 
meteorological parameters. 


Table 4. Correlation coefficients among aerosol concentrations and 
meteorological parameters 


Meteorological parameters Ion concentration 
Temperature 0.26 (0.04) 
Relative humidity -0.11 (0.38) 
Atmospheric pressure -0.13 (0.37) 
Wind speed -0.41 (0.01) 


Significant levels of the correlation coefficients are given by the values in 
brackets. The correlation coefficients among total ion concentrations and 
tempreature are statistically significant (p<0.05). Statistically significant 
correlation coefficients were found also in the case of ion concentrations and 
wind speed values. Correlation coefficients between ion concentrations and 
relative humidity and atmospheric pressure were even lower, respectively - 
0.11 and -0.13. Their significance level is higher than 0.05, showing that their 
correlations are not statistically significant. These values indicate that these 
meteorological parameters play a minor role on ion concentrations. The 
aforementioned data are based on the statistical interpretation of the 
measurement data. 


3.2. Aerosol concentrations 


The table 5 reports the average values and their standard deviations of the 
aerosol concentrations of the three modes. Major emissions occur during the 
summer season. Aerosol emissions belong primarily to fine mode. The average 
concentrations of fine aerosols during summer season were more than 1.5 
times higher than during spring and autumn seasons. 

Summer marks intensive combustion processes, construction activities, 
biomass burning and traffic (Jayaratne et al., 2010). These processes are 
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predominant aerosol sources in summer; construction activities and biomass 
burning occur also during other seasons. The presence of coarse aerosols is 
related to local wind re-suspensions processes. 


Table 5. Seasonal-averaged aerosol concentrations 


Summer |__ Spring and Autumn 
Fine (10° m”) 5.2+0.9 5.0+0.7 
Mixed (10° m°) 1.1+0.6 0.9+0.5 
Coarse (10* m”) 10.4+8.6 1.9413 


To estimate the effect of meteorological parameters on aerosol concentrations, 
their Pearson correlation coefficients were determined (table 6). 


Table 6. Correlation coefficients among aerosol concentrations and 
meteorological parameters 


| Fine Mixed Coarse | 
Temperature 0.06 (0.63) 0.14 (0.27) 0.47 (0.01) 
Relative humidity 0.80 (0.01) 0.73 (0.01) 0.22 (0.08) 
Atmospheric pressure 0.20 (0.16) 0.52 (0.01) 0.01 (0.99) 
Wind speed 0.16 (0.22) 0.28 (0.03) —0.18 (0.17) 


Relative humidity has a greater impact on aerosol concentrations, especially 
those of the fine and intermediate modes than the other meteorological 
parameters. The instrument here used measures also liquid aerosols which 
concentrations increase under the high relative humidity conditions. On the 
other hand, under high temperature conditions earth surface undergoes on low 
soil moisture values. Under these conditions, re-suspension processes are 
favored, and aerosol concentrations of coarse mode are enhanced. Statistically, 
no significant correlation was obtained in the case of the wind speed and coarse 
aerosol concentrations. This was due to the limited number of available 
measurement data. 


3.3. Binom ion-aerosol 


Here, the correlation between ion and aerosol concentrations is investigated. It 
is a well known the fact that fine aerosols (because of their major number) 
influence more on the reduction of the concentrations of cluster ions, expressed 
by eq.1 (Dhanorkar and Kamra, 1992; Reischl et al., 1996; Makela et al., 2001). 

To estimate this contribution in terms of correlation coefficients, should be 
determined these coefficients (Hoppel and Frick, 1986). The correlation 
coefficient between the concentrations of cluster ions and fine aerosol particles 
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is -0.21 (0.15)—statistically insignificant. A correlation of -0.31 (0.03) was 
found in the case of the mixed aerosols. Coarse aerosols were weakly correlated 
with cluster ions with correlation coefficient 0.23 (0.12). Thus, the presence of 
fine and intermediate aerosols at the site affects the populations of the cluster 
ions, reducing their concentrations. 

Table 7 reports correlation coefficients among ion concentrations, aerosol 
concentrations (mixed mode) and meteorological parameters. 


Table 7. Correlation coefficients among ion concentrations and principal 
influencing factors 


Mixed aerosols Temperature | Wind speed 
-0.31 (0.03) 0.26 (0.04) -0.41 (0.01) 


The data in Table 7 reports that wind activity, aerosol concentrations and the 
temperature are of great impact for the ion concentration. The concentration 
of the aerosols decreases ion concentrations due to the attachment processes 
between them. Wind activity decreases the ion concentration, due to the 
dispersion process. Meanwhile, temperature acts indirectly, increasing the radon 
concentration and consequently the concentration of cluster ions. Anyway, for 
an accurate determination of the influence of meteorological parameters on 
the ion concentrations, the mutual relationship between ion and aerosol 
concentrations as well as the meteorological parameters should be taken into 
account. In addition, a larger set of data will enable to obtain more significantly 
correlation coefficients. 

Figure 2 depicts ion concentrations and influencing factors in different seasons. 
The correlation coefficients of these two sets of data (summer and other 
seasons) are very high, 0.97. This is another important result, which shows 
that these factors undergo almost in the same direction regardless of the season. 

All units of the quantities in figure 2 are expressed based on the tables 2, 3 
and 5. 
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Fig. 2 Ion concentrations and influencing factors in different seasons. 
4. CONCLUSIONS 


The data of ion and aerosol concentrations at the Velipoja site are here 
reported. The average values of these concentrations indicate generally fresh 
air, (1000-5000 cm). The source is anthropogenic activity (traffic, combustions 
and vegetation). The source for fine aerosols is also anthropogenic activities 
(traffic, combustion and construction). Re-suspension processes mostly affect 
the coarse mode. 

In addition to anthropogenic activity, the meteorological parameters have a 
great impact on these concentrations. Temperature has a greater impact on ion 
concentration. Relative humidity has a great impact on the concentrations of 
fine and mixed aerosols. Wind activity has a great impact on the concentrations 
of the coarse aerosols and the air ions. Mixed aerosols reduce the concentrations 
of cluster ions. 

Thus, principal sources and eliminators of cluster ions on a seashore site are 
determined, within a wide range of period, based on the comparison of 
measurement data during summer, spring and autumn. 
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ABSTRACT 


A survey was carried out to investigate the employment opportunities of Albanian 
newly graduated students and the data are here reported. The data mining classification 
techniques such as NaiveBayesUpdatable, Logitboost and Adaboost1 algorithms were 
in the present investigation applied. The data collected were investigated based on: 1) 
prediction accuracy, ii) model building time and, iii) error rate. In the end, the most 
appropriate algorithm was determined. 

Keywords: Data mining, Classification Model, Graduates Employment, Prediction 


1. INTRODUCTION 


Elevated unemployment imposes significant costs on the individual and society. 
The present paper investigates the employment rate of Albanian newly graduated 
students. In Albania, unemployment rate is high and very concerning. Given 
the social-economic situation, such information would be relevant for the 
institutions involved in the area. 

The quality of education system is very important for the employment of 
graduated students. Cela and colleagues (2015) reported that number of students 
that fail to graduate is high. Consequently, finding a job is hard. Given the 
situation, categorizing the students’ profile would be very important. As this is 
anew and unexplored investigation area in Albania, difficulties in data collection 
even from a non-computerized system were encountered. 

The present investigation aims at predicting the success of the graduated 
students via data mining techniques. In the end, the data mining technique was 
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determined based on the accuracy. 
Figure 1 depicts the process here carried out. 
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Fig. 1: Knowledge discovery process 


2. MATERIALS AND METHODS 


In the present investigation, the efficiency of three data mining techniques 
such as neutral networks, the decision trees and the Naive Bayes were 
compared. The data were collected from a survey sent via various social 
networks to a group of graduated students whose demographic data were 
collected along with their education and employment experience. 


2.1 Literature review and definitions 

Data Mining is the process that uses several tools for data analysis to find 
models and connection between the data, in order to make the prediction (Witten 
and Frank, 2005). 

Regarding the academic process, data mining could be used to: 1) estimate 
students’ performance, ii) make recommendations for them, 111) estimate the 
know-how process, ili) prepare them to adopt and use different techniques in 
their job based on job market tendencies and, iv) foster intellectual pursuits in 
different areas. Consequently, the job market wouldn’t be challenging (Romeo 
and Ventura, 2007; Al-Radaideh and Nagi, 2012). 

The future of a country depends on its youth potential which source is the 
education system (Muntean et al., 2010). 

Data mining is used to estimate the employment opportunities, unemployment 
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rate, job market requirements etc. In addition, it could be used by different 
companies when hiring. The cause-consequence relationship is a means to 
address education system towards job market tendencies and development 
practices towards multiple investments (Muntean, et al. 2010). 

As classification models are constructed on the basis of historical data, data 
mining might lead to racial, social, religious discrimination etc. Consequently, 
accuracy would be relevant. 


2.2 The study process 

Zhang and colleagues (2010) said that evaluating the results of students into 
the job market involves certain criteria. In the present paper, data on education 
level have been investigated. 

The present investigation provides some information about the employment 
opportunities of the newly graduated students in Albania. Consequently, different 
methods were applied and in the end the most appropriate one was determined 
based on prediction accuracy. 


2.3. Data collection 

The data were collected from a survey sent via various social networks to a 
group of graduated students aged between 18 and 27 in Tirana and Korga. 
Their demographic data were collected along with their education level, working 
experience, age, gender, residence job selection criteria etc. This age group 
was chosen as they are expected either to attend university or foster intellectual 
pursuit. The survey aimed to give accurate predictions on the employment 
possibilities of the Albanian newly graduated students. The unemployed 
graduated students were asked on the reasons of unemployment. The survey 
was sent online via social networks such as Facebook, Google+, LinkedIn, etc. 
Once the incomplete data were removed, the total number of the newly graduated 
students involved properly in the survey was 180. 

The model of success of graduated students was created. Here, success of 
the model was evaluated based on their employment success (based on the 
result variable-employed- yes/no). The prediction model consists of 10 variables 
and the Employed class. The chosen variables are in the Table 1 reported. 
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Table 1 The Chosen Variables 


Variables Description Values 
Gender Student Gender {M,F} 
Age Student Age 18-27 
Residence Student Residence { A---Tirana, B---Korca} 
{A —PhD, B- Bachelor, 
Completed Type of studies completed or C- Master, D — High School, 
Studies that are attending E- Primary School, 
F- Attending Bachelor/Master } 
Country of Country where you have { A--Albania, B--Abroad} 
Studies studied 
Type of School | The type of school attended {A---Public, B---Private} 
Program Study Program 
ie The last study program Average 
eee (or the average of Higher Education, if is still attending) 
Foreign The number of the foreign {A---1, B---2, C ---3, D--- >3} 
Language languages you can speak 
eins Have you attended training or {A---Yes, B---No} 
Training 
supplementary courses 
Employment Your employment status { A---Employed, B---Unemployed} 


2015/Vol. XX (2) 


2.4. The Data mining methods 

In the present investigation The Naive Bayes Updateable (NBU), The 
LogitBoost (LB) and, The AdaBoost! (AB1) classifier were involved. 

The NBU classifier is part of Naive Bayes classifiers, which are based on 
Bayes’ rule and “naively” assumes independence—they are only valid to multiply 
probabilities when the events are independent (Witten and Frank, 2005). 

The LB classifier was developed by Friedman and colleagues (2000) and 
performs an additive logistic regression and generates individual models in every 
interaction. The LogitBoost algorithm uses the cross-validation to determine 
an appropriate number of iterations and finish its execution before the 
performance comes down (Witten and Frank, 2005). 

The AB1 classifier is an algorithm which is specifically designed to classify 
and which can be used also as a model in every classifying algorithm. This 
algorithm can handle weight instances, where the weight of an instance is a 
positive number. The weighted instances treatment changes the method of 
mistake calculation (Witten and Frank, 2005). 

In addition to NBU, LB and AB1, the WEKA was used for data mining. 
WEKA was developed at the University of Waikato in New Zealand, and the 
name stands for Waikato Environment for Knowledge Analysis. This package 
has been implemented in the software language Java, and today stands out as 
probably the most competent and comprehensive package with algorithms of 
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machinery learning in academic and nonprofit world (Machine Learning Group, 
2013). As a public data mining platform, WEKA collects many machine learning 
algorithms to mine data, including data pretreatment, classification, regression, 
cluster class, association rule mining and visualization on new interface. WEKA 
is very important for the academic and industrial researchers, and it could be 
used for teaching purposes. 

Figures 2a and b depict the nominal and numeric variables of the newly 
graduated students here involved. 


Selected attribute 
Name: Gender Type: Nominal 
Missing: 0 (0%) Distinct: 2 Unique: 0 (0%) 
No. Label Count Weight 
1| Female | 107 | 107.0 
2|Male 73 |73.0 


(a) Nominal Variable - Gender 


Selected attribute 
Name: Age Type: Numeric 
Missing: 0 (0%) Distinct: 9 Unique: 0 (0%) 
Statistic Value 
Minimum | 18 
Maximum |27 
Mean 23.761 
StdDev 1.732 


(b) Numeric Variable - Age 


Fig. 2 Information displayed in the Selected attributes section 


The cross validation on the classifier is very important for the evaluation of 
the classifier’s robustness (Han et al., 2001; Kantardzic, 2011). In the present 
investigation, a 3-fold cross validation consisting of a randomly split data set 
into 3 subsets of equal size was used. Two subsets were used for training. One 
subset was used for cross validating, and one for measuring the predictive 
accuracy of the final constructed network. This procedure was performed 
three times so that each subset was tested once. Test results were averaged 
over 3-fold cross validation runs. Data splitting was carried out without sampling 
stratification. The WEKA software toolkit can calculate all these performance 
metrics after running a specified k-fold cross-validation. The prediction accuracy 
of the models was compared. 


3. RESULTS AND DISCUSSIONS 
Once the chosen algorithms were executed, the results were assessed and 
analyzed and the data are here reported clearly. 
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3.1 The input variable analysis 

Graphs for each input variable in relation to the employed class are in the 
Figure 3 depicted. In these graphs, the dark part shows the level of the 
unsuccessful group (in this dataset it shows the unemployment), while the lighter 
part shows the successful group. In the present investigation 180 newly 
graduated students are involved. 100 students were from Tirana and the 
remainder from Kor¢a. The following data could be obtained from the graph 3: 
1) Tirana marking the higher employment rate, ii) the young generation today 
prefers to attend the university. 17 out of 180 students dropped out of education 
level. 75 attend bachelor classes, 66 students attend master degree classes 
and 2 are PhD students. The 2 PhD students are already employed. In addition 
to PhD students, the master degree students are employed as well, iii) the 
employment-education relationship shows that most of the students that have 
finished the high school are employed, iv) most of the students attend the 
Albanian universities (155 students). However, the students studying abroad 
have more employment opportunities, v) the number of students attending private 
universities is higher (58). In addition, they have more employment opportunities. 
Grade average merely has any impact on employment opportunities, vi) few 
students foster intellectual pursuits (training, courses etc), vi) the more foreign 
languages, the more employment opportunities and, vii) employment rate is 
high (126 graduated students). 
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Fig. 3: Graphs for each input variable in ration with the employed class 


3.2 Evaluation of the algorithms 

In the present investigation, three data mining algorithms were involved to 
evaluate the performance and the efficiency of the algorithms here involved 
and the results are in the table 2 reported. 

Table 3 reports the data obtained from the three algorithms— true positive 
(TP), true negative (TN), false positive (FP), false negative (FN), precision 
and the recall (Witten and Frank, 2005). 

TP rate is the proportion of examples which were classified as class X among 
all examples which truly have class x. In the present investigation, it represents 
the ratio between the data predicted as positive and all the data included in the 
“Employed” class. 
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TN represents the ratio between the data predicted as negative and all the 
data included in the “Employed” class. 

FP rate is the proportion of examples which were classified as class x, but 
belong to a different class, among all examples which are not of class x. In the 
present investigation it represents the ratio between the data classified by 
mistake in the “Employed” class and all the data included in all other classes, 
except the “Employed” class. 

FN is the number of examples predicted negative that are actually positive 

Precision is the proportion of positive predictions that are correct. In the 
present investigation it represents the ratio between all the examples that 
indeed belong to the “Employed” class and the total number of the examples 
classified as belonging to the “Employed” class. 

Recall is the proportion of positive labeled instances that were predicted 
as positive: 


F-measure a combined measure for precision and recall: 


2+ Precision * Recail 
F — Measuare = ———— 
Precision + Recall 


The performance of the three models was evaluated based on: i) prediction 
accuracy, ii) the time to build the model and, iii) the error rate (Table 2). 


Table 2 The evaluation of the performance algorithms 


Classifier 

Evaluation Criteria NBU LB ABI 
Time taken to build model (Sec) 0 0.03 0.01 
Correctly Classified Instances 133 132 133 
Incorrectly Classified Instances 47 48 47 
Accuracy of prediction 73.89% | 73.33% | 73.89% 
Kappa Statistic 0.3143 0.3324 0.3614 
Mean Absolute Error 0.3386 0.346 0.3576 
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Table 3 Detailed accuracy by Class 


Classifier | TP FP Precision Recall F-Measure Class 
0,833 | 0,481 | 0,802 0,833 | 0,817 Yes 
NBU 0,519 | 0,167. 0,571 0,519 | 0,544 No 
0,865 | 0,574 | 0,779 0,865 | 0,820 Yes 
LB 0,426 | 0,135 ‘| 0,575 0,426 ‘| 0,489 No 
0,865 | 0,556 | 0,784 0,865 | 0,823 Yes 
ABI 0,444. 10,135 ‘| 0,585 0.444 | 0,505 No 


The Confusion Matrices are very effective in the evaluation of the classifiers. 
The columns stand for the prediction and the lines for the current classes. The 
confusion matrix on the Table 4-6 showed that NBU, LB and AB1 algorithms 
provided relatively good results. As a result, data mining can help in the prediction 
of the youth success in the job market (employed: YES or NO). 


Table 4 NaiveBayesUpdatable Confusion Matrix 


NBU Predicted Class | 
A-Yes | B-No_ | 
Actual A -— Yes 105 21 
Class B- No 26 28 
Table 5 LogitBoost Confusion Matrix 
NBU Predicted Class | 
A-Yes B-No | 
Actual A-— Yes 109 17 
Class B- No 31 23 
Table 6 AdaBoost1 Confusion Matrix 
NBU Predicted Class 
A-Yes B-No 
Actual A — Yes 109 17 
Class B- No 30 24 


3.3 RESULTS 

In the present paper, three important Data Mining algorithms were involved 
to evaluate the preprocessing of the data and their performance. In the end, 
and the most appropriate algorithm was determined based on the accuracy. 

Results reported that: i) the LogitBoost classifier was the algorithm with the 
lowest accuracy (72%) and the longest execution time (0.03 sec), compared to 
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the other two algorithms. Consequently, it is inappropriate for application, 1i) 
the AdaBoostM1 algorithm had the same accuracy (73.88%) as the 
NaiveBayesUpdatable, but its execution time was longer (0.01 sec). The 
AdaBoost1 algorithm managed to classify accurately more instances of the 
class “Yes” (109) comparing to the NaiveBayesUpdatable algorithm, which 
managed to accurately classify only 105 instances from the 126 examples 
belonging to the class “Yes”, ii) based on the Table 2 and 3, the 
NaiveBayesUpdatable algorithm not only has a shorter execution time, but it 
offers higher accuracy, especially for the examples of the “No” class in which 
the other two algorithms failed. In addition, the average value of TP, Precision, 
Recall, the F-Measure and MCC were higher compared to the other cases, 
while the value of FP (which is not a positive indicator) is smaller compared to 
the classification obtained from the AdaBoost1 algorithm. Even though the 
algorithms have the same accuracy of 73.8889%, the data obtained from the 
output make the NaiveBayesUpdatable the most appropriate algorithm. In 
addition, it can simultaneously be correct and user-friendly for its users. 
Moreover, it could be used to: 1) assess the employment opportunities, ii) the 
unemployment rate, iii) identify the job market requirements, etc. It could be 
used by the companies when hiring. The cause- consequence relationship is a 
means to address appropriate education policies. 


4. CONCLUSIONS 


Graduates employability remains a national issue due to the increasing number 
of graduates produced by higher education institutions each year. Consequently, 
the present investigation aims at: 1) assisting education state institutions in 
encouraging students to foster their intellectual pursuits, ii) determine the factors 
effecting employment, iii) predicting whether a graduate has been employed or 
remains unemployed after graduation based on a series of classification 
experiments involving various algorithms under Bayesian and Meta methods to 
classify a graduate profile as employed or unemployed. Results showed that 
the NaiveBayesUpdatable algorithm, a variant of the Bayes algorithm, provided 
the highest accuracy -73.89%. The average accuracy of the two other algorithms 
was 73.61%. 

The cause-consequence relationship is a means to address appropriate 
education policies based on job market tendencies. 

As the data here reported are relevant for the state institutions involved in 
the area, the investigation must further by collecting a real and large data set 
from different sources such as graduate profiles from university student 
database and apply the model using such data. In addition, as the correlation 
factor between current and previous employers or the reason of choosing one 
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job to another will be included. Moreover, other classification methods could 
be applied to test the most suitable method for the structure of the graduate 
employability in Albania and give better classification accuracy. 
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ABSTRACT 


Albania is experiencing some geomorphological changes of great impact for the 
environment due to anthropogenic activity. These changes are considerably reflected 
in the tectonic zone of Pre-Adriatic Depression. Here, the area is characterized by 
erosion, landslides and movement of the Adriatic Sea shoreline. 

Anthropogenic activities have been intense in the hilly ranges running parallel to the 
sea shore which partly consist of lousy sandstone due to some catastrophic landslides 
occurred. Complicated geomorphologic phenomena could be noted in the sea shore 
which is characterised by significant sea transgression in typically accumulative areas. 
These phenomena have been investigated for a long time using satellite imagery 
involving free archives of LANDSAT [Kotor], MODIS [Balwois], SAR ERS and 
ENVISAT. In the present paper, some results recently obtained using SAR images are 
reported. 


1. REVIEW OF GEOMORPHOLOGIC PHENOMENA IN 
ALBANIAN LOWLANDS 


Half of Northern of the Western Albania— the Pre-Adriatic Depression 
tectonic zone—is characterized by lowlands and hilly ranges, bordered in the 
West by the Adriatic Sea (Geological Map 1:200000 of Albania; Aliaj et. al., 
2001; Favreto et. al., 2013; Frashéri et. al., 2014). Once a seaside, quaternary 
deposits consisting of lousy sandstone, conglomerates and clay are recent 
geologic formations are located today 10-15 km inshore. Drini River and Mati 
River flowing in North and Shkumbini, Semani and Vjosa River flowing in 
South have been of great impact for the development of lowlands. 
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In the area of economic activities occurred in the Adriatic shore areas, the 
last twenty years mark a significant increase of uncontrolled buildings’ 
construction in the hilly ranges harming the environment. Consequently, sea 
transgression has occurred and landslide intensified (Aliaj et. al., 2001; Frashéri 
et. al., 2014). 

False colour LANDSAT images in the figure 1 represent respectively the 
southern (a) and the northern (b) sections of Albanian lowlands. Images are 
obtained by the combination of NIR bands from three different years 1972- 
1986-2002 as a single RGB image, where sea regression and transgression 
areas are visible in Blue and Red colours, respectively (Colesanti et. al., 2003(a); 
Mannunta et. al., 2003; Catani et. al., 2006; Farina et. al., 2006; ESA 2000- 
2015(a); ESA 2000-2015(b); Tamburini et. al., 2011) 


(a) (b) 
Fig.1. Lowlands from NIR bands of LANDSAT. 
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The Investigation was carried out in two “layers” in the following areas: i) 
the delta of Shkumbini River and the Durrési region where the SqueeSAR™ 
software was used and multimage SAR calculations were carried out by the 
company TRE in Milano, Italy. Traditional two-image interferograms were 
calculated for both segments of lowlands as presented in figure | and, ii) Synej. 

Figure 2 depicts a typical landslide which occurred in 2007 in Synej. These 
hills geologically consist of lousy sandstone. Anthropogenic activities have been 
of great impact for their structure. The sliding mass was separated from the 
hill forming a deep “canyon” that now is used as a road. 


Fig. 2: Landslide in Synej. 


Figure 3 depicts the section of the Adriatic beach in Semani delta. Here, 
water tower now situated offshore could be noted. The source of the sea 
transgression is the sea erosion and ground subsidence. 


(a - 2007) (b - 2010) 
Fig. 3: Sea transgression in Semani area. 


Neki FRASHERI, Salvatore BUSHATI, Gudar BEQIRAJ, 
54 Mentor LAMAJ 2015/Vol. XX (2) 


Once traditional images from the field and satellites LANDSAT AND MODIS 
were applied, SAR images and interferometry techniques were used to evaluate 
vertical movements of the ground with centimetric precision (Avallone and 
Zollo 1999; Osmanoglu, 2010; UGS, 2010; Moisiu et al., 2014). The methods 
applied are in the forthcoming section reported. 


2. THE SAR METHOD 


The SAR differential interferometry involving the surface interferometry and 
scatterer-based point-wise interferometry are in the present investigation used. 
The EOLISA software was used to access the database and select pair of 
images for both southern and northern segments of lowlands. 

This method was applied based on (Ferreti et al., 2001; Allievi et al., 2003 (a) 
Alhievi et. al., 2003(b); Colesanti et al., 2003; Ubisci et al., 2005; Bottero et 
al., 2006; Albaviator, 2007; Reigber et al., 2007; Ruiz et al., 2008; Ferreti et 
al., 2011; Favretto et al., 2013; Giannico 2013). 

For surface interferograms formally suitable pairs of ASAR images were 
obtained from the ESA in the framework of a small project “Seman 14921”, in 
2013, to select pairs with the largest time difference possible. For the southern 
segment two pairs of IMS images were selected for the dates 2003/07/04: 
2007/02/23 with baseline 40m and 2003/03/21: 2007/11/30 with baseline 746m. 
For the northern segment only one pair of IMS images was selected for the 
dates 2003/03/21: 2004/11/05 with baseline 100m. In addition, ERS images 
with larger time difference were downloaded. The problems with the orbit files 
prevented us from their use. 

The calculation procedure was tested using the ASAR IMS images from 
L Aquila earthquake of 2009. Two software packages were tested, firstly RAT 
& IDIOT and finally ESA NEST. 

Two final interferograms were calculated using ESA NEST software with 
pairs of dates 2003/07/04: 2007/02/23 and 2003/03/21: 2007/02/23. The image 
of 2007/02/23 was not used because of related orbit file problem. Both wrapped 
and unwrapped interferograms were calculated. Interpretation of results is 
presented in the next section. 

Point-wise scattered interferograms were calculated with the assistance of 
the T.R.E co., Milano, Italy. Here, the software package SqueeSAR was used 
(Fig. 4). 
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Fig. 4. The investigated area via SqueeSAR. 


Two sets of ERS images were used for the period running from 1992 to 
2000: i) 31 images of upward pass and, ii) 50 images of download pass. Time 
interval between two successive images was 35 days. A total of 11.657 scattered 
points was identified. Both horizontal and vertical velocities of scattered points 
were calculated from the Line-of-Soght (LOS) relative displacements. 


3. RESULTS OF POINT INTERFEROGRAMS 


Figure 5 depicts the points selected from the area covered by SqueeSAR for 
detailed analysis. 

The most critical section of this area is located in the hilly range bordered 
with red at the centre of the image. Figure 6 depicts the numerous landslides 
identified in this section using traditional geological field work. 
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Fig. 6. Landslides identified through field work. 
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Annual average values of displacement velocities in LOS are in the figure 7 
depicted. 


Fig.7: Annual averages of velocities in LOS. 


Identified scattered points resulted mainly in urban areas along with those 
with scattered vegetation. Significant displacements are in Kavaja city, delta 
of Shkumbini River and the village of Térbufi. Based on LOS displacements 
and orbital parameters the vertical and horizontal displacement velocities were 
calculated (Fig. 8 and 9). 
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Fig. 8. Annual averages of vertical velocities 
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Fig. 9: Annual average of horizontal E-W velocities. 


The area bordered with red line has been subject to intense field work for 
identification of landslides and field results were compared with the SqueeSAR 
data. Satellite data helped identify new landslide sites in villages Sterbeg (Y — 
43 74 489, X — 45 54 987; Fig. 10), Kazie-1 (Y — 43 73 673, X — 45 55 103; 
Fig.11) and Kazie-2 (Y — 43 73 617, X—45 55 199; Fig. 12). 
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Fig. 10: Landslide in Sterbeg. 
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Fig. 11: Landslide in Kazie-1. 
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Fig. 12: Landslide in Kazie-2. 
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4. RESULTS OF SURFACE INTERFEROGRAMS 


The Figure 13 depicts the first differential interferograms calculated for the 
southern segment of the PreAdriatic Depression lowlands. The coherence is 
in the Figure 14 depicted. Few fringes which cannot be directly interpreted as 
vertical displacements of the ground could be met. For most of them the source 
is the vegetation variations, mostly in mountainous areas. Only fringes in the 
seashore in the point B and C are formally correlated with sea regression and 
transgression. The point A represents a wetland near the seashore. In Polis 
mountain range vegetation, variation is already identified via Landsat images. 


Fig. 13: Southern differential interferogram. 
! 


Fig.14. Southern interferogram coherence. 
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The unwrapped interferogram is in the figure 15 depicted. The image was 
filtered with Gaussian filter and a number of straight features correlating with 
the principal probably tectonic mountain ranges, fractures were identified. 


Fig. 15. Unwrapped filtered interferogram. 


The differential interferogram of northern segment of Pre-Adriatic Depression 
lowlands is in figure 16 depicted. 


pe weer 


Fig. 16. Northern differential interferogram 
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Correlation of fringes with mountain ranges is more visible than in the south. 
Here, their interpretation relates to vegetation variations. Nevertheless, weak 
signs of correlation with the slow but huge landslide of Rragami in shores of 
Vau i Dejés Hydropower Lake are identified (Fig. 17). 


(a) (b) 
Fig. 17: DEM (a) and interferogram (b) of Ragami landslide. 


5. CONCLUSIONS 


Field observations in the PreAdriatic Depression lowlands reported that 
residential buildings are constructed in the Semani area, in a dynamic terrain 
being developed. Consequently, the buildings are partly destroyed as result of a 
complex of factors that include movements of deltas, possible subsidence of 
the ground evolved in several decades, and disastrous landslides. 

There is a lack of scientific information regarding local ecosystem processes 
and inappropriate policies have been a continuous obstacle for suitable 
environmental solutions (Moisiu et al., 2014). 

LANDSAT images helped identify the areas where the sea shore line has 
moved mainly inland, but the images resolution is inadequate for a detailed 
analysis of the shore and landslides in hilly areas. Radar imagery ERS and 
ENVISAT have resolution similar to the LANDSAT. They help evaluate the 
surface and point-wise of vertical movement of the ground in areas undisturbed 
by vegetation variations. 

SqueeSARTTM and NEST software was used to develop interferograms. 
Multi-image pointwise interferograms built by the T.R.E. co. involving the 
SqueeSAR TM for a limited segment of lowlands allowed the identification of 
ground displacements related to landslides including unknown landslides ina 
studied area. The latter has been carefully scanned during the field works. 
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Subsidence was identified in some of the urban areas and in the delta of 
Shkumbini River in particular. Periodic use of point-wise interferograms for 
the entire territory would help to identify unknown, existing and potential 
landslides. In addition, appropriate measures for mitigation of negative impact 
could be taken. 

A minimal number of SAR images provided by the NEST software were 
used for surface differential interferograms, which provided few fringes in 
hilly and mountainous ranges along with potential fringes to be caused by 
vegetation variations. Only few weak fringes potentially due to beaches vertical 
movements and massive but slow landslide (in one case) were identified. 

The conclusion is simple — despite the complexity of terrain evolution in Albania 
including rapid urban expansion and significant vegetation variations, SAR 
interferograms proved to be a tool for analyzing the existing and potential 
landslides, as well as areas with significant vertical movements of the ground. 
Application of Remote Sensing technique involving the SAR intereferograms 
for the entire territory of Albania is a means to address an accurate management 
of catastrophic phenomena. 
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ABSTRACT 


The present paper reviews some of the remote sensing studies carried out to 
investigation environmental changes in Albania using simple methods and image 
processing software for general use and exploiting free internet archives of satellite 
imagery. Data of significant importance were obtained about the hot areas of 
environmental changes such as sea coasts experiencing sea transgression, temporal 
variations of vegetation and aerosols, lakes, landslides and regional tectonics. 


1. INTRODUCTION 


Satellite imagery is widely used to investigate and map environmental problems 
and changes, reducing at the same time the expensive field work. 

Remote sensing is a means to address a wide range of environmental, climatic, 
atmospheric and oceanography studies. In the present paper, information about 
some environmental hot spots in Albania regarding variations of water bodies 
and sea shorelines, evolution of vegetation in years due to anthropogenic activity 
is reported and global situation of aerosols in time, delineation of regional 
tectonics related with destructive earthquakes etc., and in particular lake areas 
and Adriatic shore line is evaluated. 


2. MATERIALS AND METHODS 


LANDSAT and MODIS images were obtained from the online archive of 
Geological Service of the United States. Some experiments were made using 
SAR interferometry involving ENVISAT ASAR images obtained from the 
European Space Agency online archive. The methods (Frashéri et al., 2010) 
here used were: i) combination of different bands to identify specific 
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environmental features (combining of a specific band from 2-3 dates in a false 
RGB color image, and using color differences to identify temporal change) 
and, ii) comparison of images from different dates for a single area (calculation, 
for a temporal suite of images of a specific band, of linear or quadratic trend in 
time for each pixel, obtaining gray scale images representing the spatial 
distribution of average values of the specific band, and of respective spatial 
distribution of temporal gradients of each pixel). 

Images processing involved professional software for remote sensing to cut 
images for selected areas, image processing of general use for combination of 
RGB channels, and in-house developed software for trend analysis in time 
domain for a suite of images of the same area. 


Ohrid —- Prespa Lakes area 


Fig.1: SW of Micro Prespa Lake in 2000 and 2006. 


Ohrid — Prespa Lakes represent a single surface and underground water 
body system, with water flow from Micro Prespa into Macro Prespa, and from 
the latter into Ohrid Lake in Saint Naum and Tushemisht springs. Prespa Lakes 
have experienced significant changes (Frashéri et al., 2010; Frashéri et al., 
2012; Frashéri et al., 2013). In the south-west, the Micro Prespa Lake is filled 
with sediments from Devolli River, due to artificial channels used for water 
management (Fig.1). The water surface of the Macro Prespa Lake was reduced 
several meters of height. Apparently, the source is the underground karst 
changes. In addition to field observations, satellite imagery was used to evaluate 
regional environmental changes from 1972 to present. 
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Fig.2: Ohrid-Prespa Lakes: Landsat GRNir band combination for years 1972 and 2001. 
White circles show areas with loss of free water surface and development of reeds. 


The combination of Landsat bands helped identify the dynamics of shorelines 
due to water loss in Macro Prespa Lake and sedimentation in Micro Prespa 
Lake. In the Figure 2 and 3, the combinations of Green-Red-NIr bands for 
1972 and 2001 show clearly the loss of free water surface in Micro Prespa 
Lake. Combination as RGB of NIr bands from 1972, 1986, 2002, 2010 is depicted 
in Figure 4, where the dominance of Red indicates areas with loss of free 
water surface and a slight improvement in Micro Prespa Lake in 2010 when 
the use of artificial channel ended. 


Se 
Fig.3: Micro Prespa Lake: Landsat BGNir band combination for years 1972 and 2001. 
White circles show areas with loss of free water surface and development of reeds. 
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Fig.4: Micro Prespa Lake evolution from Landsat NIr images. Left: 
Blue~1972,Green~1986,Red~2002. Right: Blue~1986, Green~2002, Red~2010. 


The trend analysis in time domain for a suite NDVI views provided by the 
Landsat images of the period 1972 — 2012 for Prespa lakes area is in Figure 5 
depicted. Results gave trend coefficients for each pixel, presented with false 
colors. The image processing software of general use was involved for the 
classification of the gradient of change in order to identify areas where vegetation 
is improving. 


Fig.5: Trend analysis of NDVI for Landsat images of 1972-2012. Left: average (Red) and 
gradient (Blue-Green). Right: gradient enhanced, positive values in Blue-Green. 


The Shkodra Lake area 

Shkodra Lake is located in northwestern Albania and is the largest lake in the 
entire Balkan Peninsula, with an area of 368 hectares. The lake is shared 
between Albania and Montenegro, of which 149 hectares and 57 km of shoreline 
fall within the Albanian territory. The lake is connected with two important 
rivers — Drini and Buna. Buna River flows from the Lake into Adriatric Sea. 
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Drini River is the source of the complicated hydrological system which 
characterizes the area because it has shifted from Shéngjini Bay to Buna River, 
due to human intervention. 
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Fig.6: Left: changes of Adriatic shoreline in the delta of Buna River, combination of 
Landsat NIr bands for years 1975 (Red) and 2009 (Blue-Green). Right: aerial view of Buna 
River delta. 


Figure 6 depicts the delineation of the erosion of delta of Buna River and 
with the accumulations in Velipoja beach obtained due to the combinations of 
Landsat NIr images from 1975 and 2009 (Frashéri et al., 2011; Frashéri et al., 
2014). 

The same method was used for the Shkodra Lake, delineating sedimentation 
areas in northern shore and near the area where Buna River flows. The source 
of the latter might be the Drini River which in winter flows into the Lake. 
Figure 7 depicts the Northern shores accumulation characterized by fast 
development of reeds beds. 


Fig.7: Left: sedimentation in Shkodra Lake — Landsat NIr bands combiantions for 1975 
(Red), 2000 (Green) and 2009 (Blue). Right: NDVI trend analysis from Landsat images, 
averages as Red, gradients as Blue-Green (visible the increase of vegetation in new sediments 
in northern shores of Lake). 
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The Drini — Buna Rivers are mentioned for periodic floods in winter. 
Combination of NIr bands from MODIS was used to delineate flooded areas 
from these rivers in 2010 and recent floods from Vjosa River in the south 
(Fig.8) 


River, Dec 2014 as Red, Jan 2015 as Blue-Green. Right: Buna River, Dec 2009 as Red, Jan 
2010 as Blue-Green 


Adriatic Sea coastal area 
The Adriatic Sea coastal areas are mostly Quaternary accumulations, part of 
Preadriatic Depression tectonic zone. 


Fig.9: Landsat NIr band combinations for southern (left) and northern (right) parts of 
Adriatic coastal zone: 1972as Red, 1986 as Green and 2002 as Blue. White linesh point 
coastal sites where sea regression or transgression is observed. 
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Adriatic coast represents about two-thirds of country’s sea coasts and is 
experiencing important changes related with regression and transgression of 
the sea, result of a combination of subsidence and erosion-accumulation that 
have strongly impacted human activities in respective areas (Frashéri et al., 
2012; Frashéri et al., 2014; Frashéri et al., 2015). Comparison of Nir bands of 
Landsat images from 1972 were a means to address the delineation of regression 
and transgression areas, in particular the movement of rivers deltas (Fig.9). 


Aerosols 

Figure 10 depicts the global temporal variation of aerosols investigated using 
the trend analysis of a suite of MODIS images preprocessed for the optical 
depth of aerosols (Pano ef al., 2004). 


Fig.10: Coverage (left) and trend analysis (right: averages as Red, gradient as Blue-Green) 
of aerosols optical depth from MODIS images for the period 2000-2010. 


Detailed analysis identified coastal areas with higher levels of aerosols (Fig.11 
left), as well as some hot spots with an increasing gradient (Fig.11 right). 


Fig.11: Trend analysis of aerosols optical depth. Left: the average (low values in Red, high 
values in Blue-Green). Right: the gradient (negative in Red, positive in Blue-Green). 


The main source is the reduction of the total of aerosols for the whole column 
of atmosphere in most of the territory of the country. This does not mean a- 
priory that distribution of aerosols near the ground surface follows the same 
trends. 
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Vegetation 

Figure 12 depicts the evolution of vegetation investigated via MODIS images 
involving the trend analysis of Normalized Difference Vegetation Index (NDVI) 
for the entire territory of the country (Frashéri et al., 2010). 


Fig.12: Trend analysis of NDVI from MODIS images (Red and NIr bands), period 2005- 
2010. Averages in Red, gradient in Blue-Green. 


As it is clearly shown, there are complicated patterns of vegetation, especially 
in mountainous areas. In general an increase of vegetation index is observed, 
related with the abandonment of mountainous areas and increase of raining in 
summers. Analysis of forestry areas is not included in this study. 


Application of Radar Interferometry 

ESA Envisat ASAR images were used for the calculation of differential 
interferograms for the North-Western Albania. Clear fringes which delineation 
was correlated with near coast hilly ranges in Shkodra area were obtained (Fig. 13). 
The source is the vertical movement of the ground and vegetation variations. 
The type of fringes delineation is related to the vegetation (Pano et al., 2014). 


JNTS Application of remote sensing for the analysis of environmental changes in Albania 77 


(b) 
Fig. 13: Environmental changes in northern Preadriatic Depression from Envisat SAR 
imagery. Left: Combination of intensity bands from 2003/03/21 and 2004/11/05. Rigt: 


DINSAR inerferogram from the same images. 


Geological Applications 
Efforts for geological applications of remote sensing were made while 
investigating the Shkodra—Peja traversal fracture (Pano et al., 2006). 


Fig.14: Combination of Landsat bands: 7 as Red, 4 as Green, 2 as Blue. Red line represent 
the Shkodra-Peja fracture. Mirdita massif (magenta area in center) distinguishes from 
limestone of Alps north of the fracture Green areas represent the vegetation coverage. 


A combination of Landsat bands 7 as Red, 4 as Green, 2 as Blue was used to 
detect for the changes, soil types and moisture (Fig.14). Tectonic zone of 
ophiolites of Mirdita is clearly separated from limestone Alps. The separation 
line coincides with the fracture as determined from other studies. 

Geophysical investigation was carried out for slope stability and landslides 
occurring in Ragami in the hilly shores of Vau Dejes artificial lake involving the 
Envisat ASAR differential interferometry (Fig.15). 
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Fig.15: Ragami landslide (black arrow ) in Vau Dejes hydroplant from Envisat SAR 
images. Left: Digital Elevation Model. Right: DINSAR interferogram. 


The landslide looks like a weak fringe in a background of vegetation changes. 
3. CONCLUSIONS 


Remote sensing satellite images are appropriate for environmental studies. A 
number of unknown environmental phenomena were identified and some known 
phenomena spatially and temporary delineated. 

The free Internet published LANDSAT and MODIS archives from USGS, 
as well as Internet archives of European Space Agency ESA accessible in 
framework of collaborative projects resulted an excellent source for regional 
studies. 

Significant data concerning the spatial and temporal environmental changes 
were obtained using image processing software for general use combining 
properly different bands as a single false color image, where different basic 
colors represent the environmental status in different conditions. 

Further complex studies blending together different spatial and temporal 
environmental parameters from satellite images, calibrated with field data would 
be necessary. 
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ABSTRACT 


The present paper aims to: i) determine the incidence of the most common NI namely 
blood stream infection (BSD, in the Paediatric Intensive Care Unit (PICU) of University 
Hospital Center (UHC) and, ii) define the risk factors associated with them and the 
appropriate prevention and control measures.We performed a prospective study on 
the incidence of BSI from May 2011 to December 2012 in a single15 bed-PICU. Six 
hundred and fifteen patients admitted to PICU, who stayed more than 48 hours were 
enrolled and monitored for BSI till discharge from ICU or death. Primary BSI was 
defined as per CDC criteria. Data of patients with BSI was compared with those without 
BSI. Outcome was measured as length of PICU stay (LOS) and survival or death. In the 
present study the use of devices (endotracheal tube [ETT], central venous catheter 
[CVC] and urinary catheter [UC]) was recorded to calculate the device associated 
infection rate (DAR), incidence density and device utilization ratio (DUR) and risk 
factors. The results reported: 1) 26 episodes of BSI, with an incidence rate 4.3 per 100 
admissions or BSI rate 5.7 per 1000 patients-days, ii) 4 clinical sepsis, 7 laboratory 
confirmed-BSI, 15 episodes of central line associated blood stream infections (CLABSI) 
among 108 patients using CVC more than 48 hours and, 111) the incidence, incidence 
density, DAR and DUR among patients using CVC was 14 %, 12 per 1000 patient days, 
30.5 per 1000 central line days, 0.39 respectively. The group age mostly affected by 
CLABSI was children under one year old. The predominant isolates were gram positive 
(n=16, 73%) of which Staphylococcus aureus was the most common, followed by gram 
negative (n=6, 27%) of which Enterobacter spp. was the most encountered. Higher 
frequency and duration of device utilization of CVC, ETT, UC, use of parenteral nutrition, 
primary surgical diagnosis, younger age, were risk factors for BSI on univariate analysis. 
The median LOS was longer in patients with BSI compared with those without (17.1 
vs.7.05, p=0.0001). However, the presence of BSI, was not associated with a raise in 


Irena KASMI, Sashenka SALLABANDA, 
82 Gentian KASMI, Eliana IBRAHIMI 2015/Vol. XX (2) 


morta- lity rate which was similar in patients with and without BSI respectively 23.1% 
vs.17.3% (p = 0.3). BSI rate in our study rate was consistent with that reported by other 

similar study. The device associated BSI was higher than international standards. This 

study has established a benchmark for future comparisons. Preventive measures such 

as improvements in central line insertion and maintenance are necessary to lower this 

baseline rate in the future comparisons. 

Keywords: nosocomial blood stream infections, children 


1. INTRODUCTION 


Intensive care units represent the largest source of infection within the hospital 
(Donotwitz et al, 1982; Donotwitz 1986). The incidence of NI in PICU varies 
from 6 to 13.7% in developed countries. In developing countries this incidence 
can be higher, up to 29.6% (Aserbergiene et al., 2009). In particular, device 
associated infections as blood stream infections (BSJ), incidence densities can 
be up to 19 times higher than in developed countries. Primary Bloodstream 
Infections (BSIs) are the main Health Care Associated Infections (HAIs) in 
the USA and the second HAIs in European PICU patients, with substantial 
effect on cost and outcome (WHO, 2011). Health care associated infections 
rate are particularly high in neonatal and pediatric ICUs. In studies in neonatal 
and PICUs in developing countries, CLABSI densities were found to range 
from 10-2 and 60 (median 18-7) episodes per 1000 catheter-days (Allengrazi et 
al., 2010). Most of the literature concerning BSI is generated from adult studies. 
Because of significant differences in age, underlying medical conditions and 
type, and distribution of pathogen these data cannot be directly extrap- olated 
to children (Cobb et al., 1993). 

Although consolidated data on the prevalence of NI in Albanian hospitals is 
lacking, it is estimated to be in the range of 10-31.6% (Faria et al., 2007; Bujari 
and Kasmi, 2009). Given the importance of the problem, we decided to study 
the incidence, risk factors, etiological organisms, and outcome of primary 
nosocomial BSI in our PICU. The results here obtained would be of great 
importance for the newly created Committee for the Infection Control of the 
hospital and Prevention already created within the hospital in implementing 
control measures against this preventable complication of intensive care. 


2. MATERIALS AND METHODS 


This was a prospective study carried out in the single PICU of the UHC, 
Tirana, Albania from May 2011 to December 2012. This PICU is a 15-bedded 
unit having an annual admission average of 450-500 patients with need for 
invasive and non-invasive treatment. All patients with an ICU stay of >48 h 
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were involved consecutively into the study. Patients, who died within the first 

24hours of admission, stayed less 48 h and those with BSI where the organism 

was the same as the one cultured prior to PICU admission were excluded 

from the study. Clinical and laboratory criteria for definition of primary BSI 

were based on the CDC criteria for infections (Garner, 1996). Criteria applied 

to patients with >48 h of PICU stay and when the organism cultured was not 

related to infection at another site. Although all consecutive PICU admissions 

during the study period that required a CVC, were prospectively studied. In 

order to esti- mate the true incidence of CLABSIs, the type of CVCs were 
included: e.g. simple (not tunneled, not antiseptic/antimicrobial) temporary 

polyurethane catheters. Insertion of temporary CVCs were done under 

anatomi- cal land mark points and maximum sterile barriers, by the attending 
physician, using the Seldinger technique. Skin antisepsis was done with povidone 

iodi- ne 10% and alcohol 70%, due to lack of the recommended 2% chlorexidine 
gluconate in Albania. No guide wire exchange was performed. Sterile gauze 

and semi-transparent patches were used as catheter maintenance due to lack 

of antimicrobial catheter dressing over insertion site. Due to the observational 

character the study did not require any deviation from routine medical care. 

The catheters stayed as long as needed if they were functioning without evidence 

of local or systemic complications. 


BSI Definitions 

An infection was defined as NI if it occurred 48 hours after admission to the 
unit or within 48 hours or 48 hours after discharge. 

Bacteraemia was defined as the biological documentation of infection, 1.e. 
the result of a positive blood culture. Laboratory-confirmed bloodstream 
infection must meet at least one of the following criteria: 1) patient with a 
recognized pathogen cultured from one or more blood cultures not related to 
another site, ii) patient with at least one of the following signs or symptoms: 
fever (>38° C), chills, or hypotension and at least one of the following: common 
skin contaminant (e.g., diphtheroids, Bacillus sp., Propionbacterium sp., 
coagulase-negative staphylococci, or micrococci) is cultured from two or more 
blood cultures drawn on separate occasions or common skin contaminant (e.g., 
diphtheroids, Bacillus sp., Propionbacterium sp., coagulase-negative 
staphylococci, o micrococci) is cultured from at least one blood culture from a 
patient with an intravascular line, and the physician institutes appropriate 
antimicrobial therapy. 

For patient <1 year of age has at least one of the following signs or symptoms: 
fever (>38° C), hypothermia (<37° C), apnea, or bradycardia and at least one 
of the mentioned above. Clinical sepsis must meet at least one of the following 
criteria: Patient has at least one of the following clinical signs or symptoms 
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with no other recognized cause, fever (>38° C), hypotension (systolic pressure 
<90 mm), or oliguria (<20 cm? /hr) and blood culture not done or no organisms 
or antigen detected in blood and no apparent infection at another site and 
physician institutes treatment for sepsis (Garner et al., 2004). 


Sampling 

Blood from a peripheral vein was drawn for culture, on clinical suspicion of 
sepsis, starting after first 48 h of ICU stay. The site of sampling was swabbed 
with povidone iodine. Once swabbed, it was allowed to dry for 2 min. Once 
dried, it was swabbed with alcohol. Two milliliters of blood was drawn and 
added 1 ml in each bottle containing 20 ml of media. Samples were delivered to 
the laboratory within an hour of collection and those collected overnight were 
stored at room temperature and delivered the next morning. 


Data collection 

Patients in ward were examined by physicians on duty who are familiar to 
the phenomenon. Patient related data, demographic details, clinical profile 
(medical, surgical, trauma), place of admission, indication for PICU admission, 
were recorded in all at the time of enrolment. 


Device utilization 

Data regarding duration of device utilization (in days) was recorded in all 
with respect to CVC, mechanical ventilation (MV) and UC. Invasive 
procedures and medications data regarding parenteral nutrition, blood 
transfusion, antibiotic use were recorded. 


Outcome This was measured as length of ICU stay (LOS) and survival or 
death. 


Statistical analysis 

The patients for analysis were divided into four age groups: i) infants (aged 1 
month to <1 year), 11) preschool children (1-5 years), iii) children (6-12 years) 
and, iv) teenager (>12 years). NI incidence (number of NI divided by number 
of patients admitted), NI incidence density (number of NI divided by number 
of patient-days), device-associated rate (number of NI divided by number of 
device-days), and device utilization ratio (number of device-days divided by 
patient-days) were calculated and compared among the age groups. Data are 
presented as mean + SD, median and percentages wherever applicable. The 
Student’ s t-test and Mann-Whitney U-test were applied for the parametric 
data. Chi-square was used for the categorical data. Data of patients with BSI 
was compared with those without BSI with respect to demographic details, 
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primary diagnosis surgical, nonsurgical, device utilization and invasive 
procedures to identify risk factors for BSI. The statistical package used in the 
study was IBM SPSS version 20. 


3. RESULTS 


A total 800 admissions were made at PICU from May 2011 to December 
2012. 615 children, 321 boys and 294 girls were in the present study involved 
as they fulfilled the inclusion criteria. The mean age of the study population 
was 19.4 months (range 0.01-168). There were 22 episodes of primary 
laboratory confirmed BSI in 22 (3.5 %) patients and 4(0.7%) diagnosed as 
clinical sepsis; the incidence being 11.2 episodes /1000 patient days. BSI was 
more common in infants less than one year of age (n=18, 81%) compared with 
children older than two years of age (n=4, 19%) (Table1).One hundred and 
eight patients required CVC insertion. Among them, 15 episodes of CLABSI 
were confirmed. The mean rate of CLABSI was 30.5 per 1000 central-catheter- 
days with a catheter utilization rate of 39% (Table 2). Of the 15 CLABSIs, 12 
isolates were gram positive bacteria, S.aureus (n=10,45%), CoNS (n=2, 9%),3 
isolates were gram negative, Enterobacter spp.(n=1,4.5%), Klebsiella pneumonia 
(n=1,4.5%), no identified gram negative (n=1,4.5%). Among the 7 laboratory 
confirmed BSI episodes 4 isolates were gram positive, CONS (n=3,13.5%), 
Staphylococcus MRSA(n=1,4.5%), 3 isolates were gram negative, Enterobacter 
sp.(n=2, 9%), Escherichia coli (n=1, 4.5%). 


Table 1 Characteristics of 615 patients in the PICU with and without BSI 
during the study period 


Characteristics Total (N =615) | With BSI | NO BSI=589 | P value 
(N=26) 
Age (months) Mean (min | 19.4(0.01;168) 19.9(0.01;144) 19.3(0.01;168) 0.92 
max) 


Sex n (%) 

Male 321(52.2) 13(50) 308 (52.3) 
female 294(47.8) 13(50) 281(47.7) 
Imo- <1 year 377(61.3) 20(80) 357(60.7) 
2-5 years 185(30.1) 4(15.3) 181(30.7) 
5-12 years 35(5.7) 0(0) 35(5.9) 
>12 years 18(2.9) 2(7.7) 16(2.7) 


Primary diagnosis at 
admission, n (%) 


Medical 437(71) 10(38.5) 427(72.5) 
Surgical 178(29) 16(61.5) 162(27.5) 0.001 
Emergency department 450(73.1) 10(38.5) 440(74.7) 


Others wards 165(26.9 16(61.5) 149(23.7) 
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Device utilization, CVC n | 108(17.5) 15(68.2) 93(15.7) 
(%) 108(17.5) 15(68.2) 93(15.7) 
Type of catheter,CVC n 0(0) 0(0) 0(0) 
Non tunnelled simple 
Permanent 

111(18) 20(76.9) 91(81.9) 
Parenteral nutrition, n(%) 208(33.8) 24(92) 184(88.5) 
Blood transfusion, n(%) 345(56) 19(73) 326(55.3) 
ETT, n(%) 339(55) 22(84.6) 317(53.8) 
CU, n(%) 

108(17.5) 6(23) 102(17.3) 
Outcome Death n, (%) 0.8(+ 2.3) 5.5(+4.9) 0.59(+1.96) 
Line duration mean (SD) TAT(+ 6.1) 17.1 (+8.6) 7.5(+5.68) 
PICU stay(days)mean (SD) 


0.001 


0.001 
0.001 
0.005 
0.001 


0.43 
0.001 
0.001 


Table 2 CLABSI incidence rates and device utilization ratio in the age groups 


Central 
Bloodstream nosocomial infections line 
Age groups Number of NIs(n) Incidence Device 
Incidence Device Density associated | utilization 
rate Ratio 

1 month - <1 year 13.0 18.5 15 40.6 0.37 
1-5 years 2.0 5.3 6 17.5 0.37 
6-12 years 0.0 0.0 0.0 0.0 0.38 
>12 years 0.0 0.0 0.0 0.0 0.81 
All 15.0 14.0 12 30.5 0.39 


Further analyses of the pathogens showed 16(73%) to be Gram-positive, 
with Staphylococcus aureus as the leading causative pathogen (N=11, 50%) 
followed by Coagulase-negative staphylococci (n=5, 22.8%). Among Gram- 
negative Enterococcus spp. were the leading causative pathogens, responsible 
for 3 (13.5%) of the 6 (27%) cases in the Gram-negative group (Table 3). 


Table 3 Etiological organisms for BSI (n=22) 


Staphylococcus aureus (S.aureus) 10(45.5) 
Staphylococcus aureus (MRSA) 1(4.5) 
Staphylococcus epidermidis(CONS) 5(22.8) 
Enterobacter spp. 3(13.7) 
Escherichia coli (E.coli) 1(4.5) 
Klebsiellapneumoniae (Kleb) 1(4.5) 
Others gram negative 1(4.5) 


Risk factors associated with BSI included younger age, primary diagnosis 
(surgical, medical; P<0.001), parenteral nutrition (P<0.001), the use of ETT 
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(P=0.005), the use of CU (P<0.001). Variables significant on univariate testing 
(P <0 .05) were entered in multiple logistic regression analysis. Patients who 
had parenteral nutrition infused through the line and blood product transfused 
were more likely to develop BSI (adjusted OR respectively: [OR= 7.9; 95% 
CI: 2.6-19.6; p=0.031]; [OR= 10.9; 95% CI: 2.04-55.5; p=0.005]). Univariate 
analysis of the outcome of BSI revealed that patients with BSI had median 
LOS longer compared with those without (17.1 vs. 7.5 days, p=0.0001). The 
crude mortality rate associated with BSI was 23.1% compared with 17.3% in 
patients without BSI (p = 0.3). 


4. DISCUSSIONS 


Carried out for the first time in the PICU of UHC, Tirana, Albania, this 
prospective study aimed at evaluating the rate, risk factors and outcomes of 
BSI among pediatric patients using the NI CDC definitions (Garner et al., 
1996). In addition, it is an appropriate means to address for future studies in 
support of the newly created committee of infection control and prevention 
within the hospital and allow for comparison of international data. The mean 
rate of CLABSI in our PICU (30.5 per 1000 central-catheter-days) was very 
high compared with the mean rate of CLABSI reported by systematic review 
conducted by WHO which varies from 2.1-4.3 in developed countries to 10.2- 
60 episodes per 1000 catheter days in developing countries (Pessoa-Silva et 
al., 2004; Allengrazi et al., 2010; WHO 2011). Richard et al., (1999) examined 
BSI in 61 PICUs in the USA and found it to be the most common nosocomial 
infection with a rate ranging from 0.0 to 20/1000 catheter-days. Stover et al., 
(2001) reported PICU BSI rates from 23 US hospitals with wide variations in 
the rates ranging from 0.0 to 18.5/1000 catheter-days. This variability depends 
on unit related parameters such as size and settings, patient related parameters 
such as severity and type of illness, site and conditions under which the catheter 
was placed. This data are missing in this study. It is difficult to explain the 
reason for our high rate as there are no previous data with which to compare. 
The source for the high rate is partly the features of our PICU’s population 
and the variability in critical care practices and local infection control. Firstly, 
the majority of our patients were young infants (61%). It is well known that 
CLABSI is more common in young infants (Chodekar et al., 2011). It is well 
demonstrated in more than one study the high rate of CLABSI in NICU patients 
(NNIS, 2001; 2003). In the USA, PICUs have more than 50% of postoperative 
cardiothoracic patients who are quite stable, with less illness severity and shorter 
length of stay, opposite to our population with postoperative patients being 29%, 
high medical emergency patient being71%. O’Grady et al., (2002) reported 
that trauma patients documented the highest CLABSI rate, and non-elective 
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admission for medical management was found as independent predictor of 
CLABSIs in the study of Costello et al., (2009) supporting our explanation. 
Secondly, the lack in nursing staff being of a ratio nurse: patient 1.5-2:1during 
the day and 3-3.5: 1 during night shift may be another reason as is well known 
that understaffing in ICU is an independent risk factor for NI in general and 
specifically for BSI which is highly dependent on hygienic behavior (Robert et 
al., 2000; Schwab et al., 2012). Thirdly, the lack of consistency in application 
of written guidelines for the access and care of CVCs in our PICU, the use of 
antimicrobial-impregnated catheters, the implementation of annual hospital wide 
hand washing campaigns, the creation of physical barriers between patient 
beds, the use of 2% chlorexidin skin disinfectant for all central lines for all 
patients of all ages have, may have contributed to the high rate. Finally, although 
it was intended to use the CDC NNIS definition for CLABSI, this was not 
strictly applied. In our study, we did not require evidence that the catheter was 
the focus of infection. It is therefore possible that some of the episodes attributed 
to the presence of CVC, could have originated from the gastrointestinal tract 
or elsewhere. 

The distribution of causative micro-organisms for BSI in our PICU is similar 
to that reported in the NNIS data (NNIS, 2001; 2003). Gram-positive organisms, 
were predominant (73%) in our population, with S.aureus (27%) being the 
leading pathogen, whereas Gram-negative organisms were less encountered 
(27%) of which Enterobacter spp. was the most encountered. 

In multiple logistic regression analysis parenteral nutrition infusion and blood 
product transfusion through the line were found to be a risk factor for CLABSI 
in our PICU as supported by similar study (Almuneef et. al., 2006; Wylie ez. 
al., 2008). This is supported by Singh-Naz and colleagues (1996) who noted 
that children receiving parenteral nutrition and antimicrobial therapy were 
associated with the higher risk of developing nosocomial infection. It has been 
emphasized in our PICU that infusion of fluid emulsion should be completed 
within 12 h, and that of blood or other blood products within 3 h. Prior antibiotic 
therapy was not studied as a risk factor for BSI in our study, during which 
approximately 98% of the enrolled patients and 100% of them with CVC were 
taking antibiotics prior to its development, although it has been mentioned in 
studies that receiving antibiotic for more than 10 days could be a risk for 
developing nosocomial infections including lower respiratory tract, BSI, and 
urinary tract (Andrea, 2005). The severity of illness of our patients was not 
reported since the PRISM III score, which is now in use, was partially in use 
during the study period. A number of knowing risk factor for BSI as the place 
of insertion, the anatomical site, the catheter insertion necessity, factors that 
should be taken in consideration for future studies. 

The mean line duration and the PICU length of stay were longer in the BSI 
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cases. The duration of catheterization is reported as independent infection risk 
factor in some studies, a fact that was confirmed by our results, showing longer 
duration of catheterization in patients with catheters who developed CLABSI 
(Wylie et al., 2008; Rey et al., 2011). However, the presence of BSI, under the 
conditions of our study, was not associated with the mortality rate which was 
identical in patients with and without BSI (23.1% vs.17.3% (p = 0.3), in 
agreement to the majority of published data (Smith 2008; O’ Grady et al., 2011). 

The major limitation of our study relates to a relatively small study center, 
which poses barriers in the generalization of our results. Although, there are 
strong points as well because international definitions we used enabled us to 
compare our data to international standards and allow us to set a benchmarking 
for future comparisons. 


5. CONCLUSIONS 


BSI especially CLABSI rate in our study is higher than the international 
standards. We set the baseline for CLABSI rate, causative microorganisms 
and risk factors as indicative infection control measures to be taken. 
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ABSTRACT 


The present paper aims to evaluate the clinical outcome of reirradiated patients with 
bone metastasis in terms of pain improvement, local control, and development of 
myelopathy. 20 patients with bone metastasis who complained pain have been 
reirradiated from 2012-2014. The median interval between the two irradiations was 17 
months (range 10-38 months). Karnofsky performance score at the time of retreatment 
ranged from 50% to 80%, and median age was 46 years (range 26-74). The first radiation 
dose varied from 30 Gy to 50 Gy with daily fraction doses of 1.8-3Gy. The total dose of 
retreatment ranged from 8-30 Gy with daily fraction doses 3-8 Gy, and the cumulative 
dose ranged from 38 Gy to 70 Gy (median 53 Gy). All patients were followed up 
prospectively and all of them tolerated well reirradiation. The median follow up after 
the last treatment was 12 months (range 8-26 months). Here, 12 patients out of 20 had 
complete pain relief and the remainder had partial pain relief. At the last follow up only 
one patient was still alive and the others were dead from systemic disease progression. 
We didn’t detect any serious acute side effect and no patient showed treatment- 
induced neurologic abnormalities. The results reported that reirradiation of spinal 
cord within the dose range we used is a feasible treatment that provides clinical benefits 
and improvements of quality of life. 

Key words: reitradiation, spinal cord, myelopathy 


1. INTRODUCTION 


Nowadays improvements of cancer therapy have led to better survival rates 
and an increased incidence of recurrences in site of previous conventional 
radiotherapy. Many of them need to be reirradated as the only treatment option 
and are not candidates for a surgical approach because of the limited indications 
(a poor performance status, co-morbidity, presence of other sites of metastases 
or localization of the metastasis in the skeleton). 

However, reirradiation still remains a challenge for the radiation oncologist 
as a result of limited evidence. The right interval between two irradiations, 
fractionation and total dose remain unclear. 
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As spinal cord is a dose-limiting organ and radiation can induce damage at 
even a small volume of it and impair normal function, we should be very careful 
not to lead to myelopathy, which is irreversible and can radically reduce the 
quality of life of the treated patient. Radiotherapy is a well-accepted treatment 
modality for localized pain with a 60-80% likelihood of overall pain relief (Wu 
et al., 2003), however higher RT doses, larger doses per fraction, and previous 
exposure to radiation could be associated with a higher probability of developing 
radiation induced myelopathy (Makbeth er al., 1996). Experimental data report 
that the total dose of the first and second radiotherapy, interval to reirradiation, 
length of irradiated spinal cord, and age of the treated animals influence the 
risk of radiation induced myelopathy (Rades et al., 2008). 

The present paper investigates the feasibility of spinal cord reirradiation, its 
effect on clinical outcome and improvement of quality of life. 


2. MATERIAL AND METHODS 


In the present investigation 20 patients with bone metastasis have been 
reirradiated from January 2012 until December 2014. 12 patients out of 20 
were men and the reminder were women with a median age of 54 years (range 
26-74 years). Except one all other patients have only bone metastasis at the 
time of retreatment. They have received multimodality treatment before and 
also have used corticotherapy. Three patients were reirradiated to the cervical 
spine, twelve to the thoracic spine and five to the lumbar spine. Patients were 
treated for different primary tumors (see tab.1), but the reason of the second 
irradiation was pain progression and refractory to analgesics. 

The total dose of the initial radiation therapy ranged from 30 to 60 Gy, with 
single doses 1.8—3Gy. All patients had initial multifractionated radiotherapy 
treatment and sixteen of them were treated with a Cobalt 60 machine and four 
with a Linear accelerator using a three-dimensional conformal irradiation 
technique. The median interval between two irradiations was 17 months (range 
10-38 months) and the total dose of reirradiation ranged from 8 to 30 Gy, with 
daily fraction of 3—8Gy. Ten patients received single dose radiotherapy and 10 
others multifractionated treatment (6 patients 5x4Gy and four patients 10x3Gy). 
Twelve patients were retreated by a Cobalt 60 machine and eight patients by a 
Linear accelerator using a three-dimensional conformal irradiation technique. 
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Table 1. Patients’ characteristics 


Characteristic Number of patients 

Age 

SOyears 7 

5Oyears 13 

Gender 

Female 12 

Male 8 

Type of tumor 

Breast cancer 6 

Prostate cancer 4 

Lung cancer i) 

Rectal cancer 3 

Renal cancer 2 

Rectal cancer 2 

Lymphoma 1 

Visceral metastasis at the time of retreatment 

Yes 

No 1 

19 

Vertebral position 

Cervical 3 

Thoracal 12 

Lumbar 5 

Time between two RT treatments 

<12months 6 
12months 14 

Radiation schedule 

Single dose 10 

Multiple fractions 10 

3. RESULTS 


Reirradiation was tolerated well by all patients. No acute side effects = 
grade 2 was seen. One patient developed dysphagia G1 and another diarrhea 
Gl. 

After the second month of retreatment twelve patients have complete pain 
relief and eight others have only partial pain relief. The evaluation of pain was 
done based on patient self assessments of pain using a numerical rating score 
from 0 (no pain) to 10 (worst imaginable pain). Also every change in pain 
medication during follow-up was recorded. 

Afterwards we followed up the patients every three months. In the next 6 
months, no one of twelve patients complain pain as a major symptom. The 


96 Elvisa KOZMA, Agim SALLAKU 2015/Vol. XX (2) 


median follow up after the last treatment was 12 months (range 8-26 months). 
At the last follow up all patients except one who is still alive complained pain, 
but the cause of death in all patients originated from progressive disease. 

No patient showed any neurologic abnormalities affecting motor or sensory 
function. They were all able to walk and were continent for urine and stool 
until the time of death. 


4. DISCUSSION 


Reirradiation of spinal cord is still under investigation. Radiation oncologists 
are very reserved about it as a result of existing concerns that reirradiation 
may induce myelopathy and severe neurologic dysfunction. Symptoms may 
develop months until years after reirradiation (Ryu et al., 2000). Nevertheless 
actually exists prospective and retrospective data that supports retreatment of 
painful bone metastasis following initial irradiation (Grosu et al., 2002, Maranzano 
et al., 2011). 

There are three clinical situations when reirradiation can be used: in case of 
no response, after partial response from first irradiation and hope to achieve 
further pain relief, and in case of recurrence of pain during follow up after 
complete or partial response with initial radiotherapy. 

In our case all patients were reirradiated for progressive pain after partial or 
complete response with initial radiotherapy. 

Although, literature shows a 12-months local-recurrence rate of 39% with 
single dose radiotherapy and of 23% for patients treated with fractionated 
courses (Rades et al., 2009) many radiation oncologist use single doses for 
reirradiation to avoid side effects. In the present investigation single doses 
were used in 10 patients. The dose of second irradiation was also conditioned 
by patient’s conditions and their life expectancy. Patients with a poor expected 
survival can benefit from a shorter course of radiation because a shorter 
overall treatment time means less daily trips to the radiation oncology department 
and positioning on the treatment couch and as a result less efforts for the 
patient and less costs for the treatment (Van den Hout et al., 2003). 

During reirradiation are used different daily fractionations and different total 
doses which make difficult the comparison between their effects in spinal cord. 
We should convert the given physical dose to biological equivalent dose (BED), 
in order to make possible the comparison between them. 

Reirradiation has a higher risk of causing myelopathy because the cumulative 
BED to spinal cord is higher (Maranzano et al., 2001, Wong et al., 1994). 

Nieder and colleagues (2006) reported that clinical data on spinal cord 
reirradiation allowed the development of a risk score that distinguish patients 
who may develop radiation induced myelopathy on the basis of three variables 
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which are cumulative BED, BED of each single radiotherapy course, and 
interval between the courses. The risk of radiation induced myelopathy appeared 
small after cumulative. 

BED § 135.5 Gy,, when the interval from the first RT and reirradiation is 
not shorter than 6 months, and the BED of each course was © 98 Gy, (Nieder 
et al., 2006) 

In our investigation we prescribed radiotherapy regimens which cummulative 
BED did not exceed 150 Gy,, and BED of each course was £ 98 Gy, to the 
tumor to the tumor and < 75 Gy, to the spinal cord. All our patients have an 
interval between the courses more than 6 months. So, according to Nieder’s 
risk categories the majority of our reirradiated patients had an estimated low 
risk of radiation induced myelopathy and indeed no one of them developed 
myelopathy. 

Our data supports the retreatment of patients with pain from bone metastasis 
following initial irradiation, especially after initial response. The three types of 
fractionation we used (single dose 8Gy, 5x4Gy and 10x3Gy) were safe and 
effective for pain relief, but patients who were treated with multifractionated 
radiotherapy showed more complete pain relief than patients treated with single 
dose reirradiation (ten from twelve patients). Nevertheless, we should take 
into account that we treated a small number of patients from which clear 
conclusions cannot be drawn. 


5. CONCLUSION 


Reirradiation of spinal cord using the fractionations and total doses we used 
is a safe treatment with very good palliative effects and no risk of complications 
for the rest of life. 

As the number of patients here in the present investigation involved was low, 
evaluating the fraction size or the total dose with optimal clinical effects and 
minimal side effects was impossible. 
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ABSTRACT 


The present paper provides information about a hybrid optimization algorithm which 
is a combination of Bat Algorithm (BA) and Hill Climbing (HC). This algorithm 
incorporates the diversification ability of BA with intensification ability of HC. The 
hybrid optimization algorithm is verified by testing on a large set of numerical test 
functions and results obtained were compared with the results obtained by bat algorithm 
(BA) and three state-of-the-art BA variants. The results reported that the performance 
of the proposed method is better than the bat algorithm. In addition, it is comparable to 
three state-of-the-art modified BA algorithms in terms of the quality of final solution 
and its convergence rates for high dimensional problems. 

Keywords: Bio-inspired algorithm, Numerical optimization, Bat algorithm, Hybrid bat 
algorithm. 


1. INTRODUCTION 


Many nature-inspired algorithms have been invented over the past a few 
decades and applied with success to many numerical and combinatorial 
optimization problems. The aim of the optimization algorithms is to find a set of 
independent variables which are maximizes or minimizes the value of dependent 
variables. Ant Colony Optimization (Dorigo et al., 1996), Particle Swarm 
Optimization (Kennedy and Eberhart 1995), Artificial Bees Algorithm 
(Karaboga 2010), Bat Algorithm (Yang 2010), and Artificial Immune System 
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(Kim et al., 2007) with increasing popularity are examples of nature inspired 
algorithms. Since the real-world optimization problems are getting more 
complicated, higher dimensioned and more dynamic, it follows that meta- 
heuristics algorithms’ popularity will continue to increase, as the finding the 
optimal solution becomes very difficult within a reasonable amount of time. 

Bat Algorithm (BA) is one of the nature-inspired algorithms recently presented 
and successfully applied to solve numerous optimization problems in diverse 
fields. It is a population-based search algorithm that was based on the 
echolocation behavior of bats (Yang 2010). BA, like Yang’s previous algorithms; 
Cuckoo Search (Yang and Deb 2009) and Firefly (Yang 2009), combines the 
advantages of existing algorithms, especially Particle Swarm Optimization and 
Harmony search. Many researches on BA have been carried out recently, due 
to its simplicity, convergence speed and population feature such as: 1) multi- 
objective bat algorithm (MOBA) (Yang 2011), ii) self -adaptive modified bat 
algorithm (SAMBA) (Baziar et al., 2013), iii) local memory search bat algorithm 
(LMSBA) (Yuanbin et al., 2013) and, iv) adaptive bat algorithm (ABA) (Wang 
et al., 2013). BA shows good results when dealing with lower-dimensional 
optimization problems, but its performance decreases for high-dimensional 
problems as it becomes more exploitative very fast and may stack in a local 
optimum. 

Hill Climbing (HC) is related to gradient ascent. Although good for unimodal 
problems, HC doesn’t guarantee the best possible solution for multi-modal 
problems as it is easily trapped in a local optimum when many local optima 
could be met. Many researchers (Renders and Bersini 1994; Lim et al., 2006; 
Yin et al., 2006; Yildiz 2009) prefer hybrid algorithms with hill climbing to increase 
the intensification of the search due to its simplicity and its high reliability for 
finding local optimum. 

This paper proposes a hybrid bat algorithm with hill climbing to increase the 
local search of BA. The hybrid is named as BAHC. The performance of the 
BAHC is extensively evaluated on eleven benchmark functions from Yao et 
al., (1999) and compared with standard Bat Algorithm and three state-of-the- 
art modified BA algorithms. 

The remainder of the paper provides information about the Bat Algorithm, 
its variants used to compare the performance of BAHC, HC, the hybrid BAHC 
and experimental results and comparison with other algorithms on optimization 
test functions. At the end some conclusions are drawn. 


2. BAT ALGORITHM (BA) 

Bats have unique ability to detect insects and avoid obstacles by using a high 
frequency sound based system which is called echolocation. Bats emit sound 
waves and listen to the returning echoes and they can draw 3D picture of their 
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surrounded environment in their brain even in pitch black dark. Bats can tell 
the shape, size, texture of a tiny prey, in which direction the prey is heading and 
even the speed of the prey by using the delay time and loudness of the response 
and the time difference between the ears. Bats have also the ability to change 
the way they emit the sound pulses (Richardson 2008; Gandomi et al., 2013; 
Topal and Altun 2014). 

The algorithm starts with the n bats being placed randomly in the search 
space. Velocity v,, frequency f, pulse rates r,, and loudness A, are initialized 
for each bat at the beginning. Bat’s position and velocity are updated by using 
(1, 2, and 3). 


Si Din Voge dae (1) 
Vi = vee + ( =a _ xf, : (2) 
a = xtt + Vv; (3) 


where B ¢ [0,1]. Then the algorithm evaluates the finesses (solutions) and 
chooses the current global best position x" at time step ¢. After these updates, 
if the bat’s pulse rate is low (it means bat is far away from the prey), with a 
high probability it will fly near the current best bat and make a random short fly 
there. If its pulse rate is high then it should be near prey and it will again with 
a high probability just make a random fly around its current position 

Kt SAY, (4) 

where €e[-1, 1] is arandom number, while A‘=< A; > is At= 1/N Das A} 
the average loudness of all bats at this time step, where N is the number of 
bats. After this fly if the bat’s position is better than the current global best and 
its loudness is loud enough to be greater than a random number, the bat will fly 
to this position and current global best will be updated with the new one. The 
bat’s pulse rate r, will be increased and loudness A, will be decreased. Then 
again finesses will be evaluated and the current global best x” will be found. 


3. MODIFIED VERSIONS OF BAT ALGORITHMS 


Here, the theory of three state-of-the-art modified versions of BA algorithms: 
Chaotic Local Search-based Bat Algorithm (CLSBA), Adaptive Bat Algorithm 
(ABA), and Novel Adaptive Bat Algorithm (NABA) is briefly reported. 


3.1 Chaotic Local Search-based Bat Algorithm (CLSBA) 
CLSBA (Chandrasekar, 2013) is a combination of the standard BA with 
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chaotic sequences generated by the logistic map. The use of chaos makes the 
frequency adaptive and more random in nature to balance the trade-off between 
exploration and exploitation. The frequency of pulse emission is modified as 
follows in CLSBA: 


f+ =n *f@ * 0 -fO) (5) 


While the step size of random walk parameters € (Eq. 4) varies within the 
range [-1, 1] in BA, it varies within the range [-SF(t), SF(t)] in CLSBA. SF is 
the scaling factor which changes dynamically according to the Rechenberg’s 
mutation rule (Rechenberg 1994). If the value of SF(t) is low, exploitation 
increases. If value of SF(t) is high, exploration accelerates. 


3.2 Adaptive Bat Algorithm (ABA) 

Wang and colleagues (2013) introduced an improved bat algorithm to solve 
BA’s premature convergence problem. In ABA, Wang and colleagues improved 
the random walk equation, but they also modified the frequency and the velocity 
equations in BA where each bat uses the same frequency increment for the 
velocity which makes bat’s flight behavior lack of flexibility. In ABA, Wang 
proposed a new method to let each bat dynamic and adaptively adjust its flight 
speed and its flight direction. ABA is targeting to increase the speed of the bat 
which is farther from the prey. The farther the distance between the bat and 
its prey (global best solution), the faster the speed flying to its prey. 

Secondly, Wang improved the random walk by combining it with shrinking 
search. Bats which are far away from the current best, make random fly near 
to the current best. The step size of random fly within the range [-1, 1] in BA 
while it shrinks in ABA as the iteration proceeds. ABA targets to increase the 
intensification of the search by shrinking search method. 


3.3 Novel Adaptive Bat Algorithm (NABA) 

NABA (Kabir et al., 2014) incorporates two techniques within BA, which 
include the Rechenberg’s 1/5 mutation rule (Rechenberg, 1994) and the 
Gaussian/Normal probability distribution to produce mutation step size. 

NABA offers new equation to generate new solution rather than doing just 
random walk. It modifies the random walk equation (Eq.4) in line 12 of original 
Bat algorithm. NABA controls the random walk step size by the variance of 
Gaussian/Normal distribution. The modified equation is as follows: 


Xow = Xoiq + €A'N(0,0) (6) 


new 


where © is the standard deviation. 
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Kabir used the Rechenberg’s 1/5 mutation rule (Rechenberg 1994) to 
adaptively change the random walk step size and pulse rate to control the 
exploration and exploitation. If the success-rate is less than 1/5, NABA is 
exploring too much and it intends to move the bats near to the best solution and 
decreases the step size. If the success-rate is more than 1/5, NABA is exploiting 
local optima too much and it increases the step size. That will help bats to 
reach the local optima faster, but the accuracy will decrease because of the 
long step size. 


4. HILL CLIMBING ALGORITHM (HC) 


Hill Climbing belongs to the family of local search and it is related to gradient 
ascent (Russell and et al., 1995). Here, we start the search from a random 
position and iteratively test new candidate solutions around the current solution. 
If there is a better solution among the candidate solutions, adopt the new one. 
This technique improves the current value of the objective function until we 
reach a local optimum. When no local move could further improve, the search 
would have reached a local optimum of the objective function (Selman and 
Gomes 2006). The pseudo code of HC technique is shown in Algorithm 2. 


Algorithm 2 Hill Climbing 

S <= some initial candidate solution 

repeat 

R= Tweak(Copy(S)) 

If Quality (R) > Quality(S) then 

S< R 

until S is the ideal solution or we have run out of time 
return § 
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Fig. 1: The bats’ positions before and after Hill Climbing on one dimensional Rastrigin’s 
function. 


5. PROPOSED HYBRID OPTIMIZATION ALGORITHM (BAHC) 


Here, we propose a new Bat Algorithm (BA) which is hybridization of the 
Bat Algorithm and Hill Climbing, called BAHC. BA suffers from the premature 
convergence problem for high dimensional problems because the algorithm 
tends to converge very fast. 

In BAHC, HC technique is used in two parts of BA. At the beginning of BA, 
the population is randomly initialized. Once initialized, HC is used to move all 
the bats to their local best. So, many local optimum traps will be detected at the 
beginning of the search. 

Figurel depicts bats movement on one dimensional Rastrigin function when 
Hill Climbing algorithm is used. The y-axis shows the fitness values of Rastrigin 
function, x-axis shows the independent variables, and small circles on the function 
represent the bats solutions. 

As it could be noted, at the beginning of the algorithm, Hill Climbing helps 
bats to move to the bottom of the basin (the algorithm is used to find the minimum 
value/fitness of Rastigin function) with a subsequent start of the bat algorithm. 

In addition, whenever there is an improvement on bat’s solution, Hill Climbing 
could enhance the improvement. So, bats will keep moving to their local optimum 
during the search with highest probability of locating the global optimum. The 
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BAHC pseudo-code is illustrated in Algorithm 1. The differences from BA are 
shown in boldface font. 


Algorithm 1 Hybrid of Bat Algorithm with Hill Climbing - BAHC pseudo code. 
Objective function f(x), x=(X,, ... , X, y 
Initialize the bat population x, Vv, ( i=1,2,...,n) 
Modify the solutions by using Hill Climbing 
Define pulse frequency f. at x, 
Initialize pulse rates r, and the loudness A, 
while(t<Max number of iterations) 
Generate new solutions by adjusting frequency, 
and updating velocities and locations/solutions [equations (1) to (3)] 
if(rand > r,) 
. Select a solution among the best solutions 
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. Generate a local solution around the selected best solution 
. endif 

. Generate a new solution by flying randomly 

14. if(rand < A, & f(x,) < f(x’) 

15. Modify the solution by using Hill Climbing 

16. Accept the new solutions 

17. Increase r,and reduce A, 

18. end if 

19. Rank the bats and find the current best x* 

20. end while 


21. Post process results and visualization 
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As aresult, BAHC differs from the original Bat Algorithm in line 3 and 15, 
where solution is improved using Hill Climbing technique. 


6. NUMERICAL EXPERIMENTS AND RESULTS 


6.1 Simulation 

In the present investigation 11 test functions were carried out for a 
comprehensive comparison (Yao et al., 1999). The benchmark functions f, - f, 
are unimodal functions. The functions f,- f,, are multimodal functions where 
the number of local minima increases exponentially with the problem dimension 
(Yao et al., 1999) that makes these functions are very hard to optimize. The 
test functions are described in Table 1. 

BAHC was compared with original BA, CLSB, ABA and NABA to 
demonstrate the effectiveness of BAHC. In our experimental studies, parameter 
settings of the algorithms are the same as in their original papers. The population 
size of the algorithms has been kept as 50 regardless of the dimension size of 
the problem. 
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Table 1. Descriptions of the test functions 


No Name Formula Global Search Domain 
Cc aeaaes 
r-eue | )= 


f, | sphere ode | Model } -t0ss,s100 | <x,<100 


Po %; 


Schwefel's 
Problem 2.22 


i Schwefel's ys 
Ss Problem 1.2 f,09 = DP x Xj y Ea like 
f Schwefel's (x) = max {|x}, 1<1< 30} P i6iie i 
oat Problem 2.21 Lie ue Sx,< 


Generalized 
f, Rosenbrock's Sfx) =X?_ 1000, —x7F + (x; - TI 0 -30 Sx, < 30 
Function 


Sfx) = L?,_, be] + Tbe — -10<x,<10 


fe Step Function FAX) = L?,_,( Le, + 0.5] P -100 < x, <100 


: Quartic Function Rat ag 
i net (pe EP, t+ ¢ 1,28 <x,< 12 
f i.e. Noise Sx) =X? _, x;1 + random [0,1) 1.28 <x, < 1.28 


Generalized 
ti Schewefel's SAX) = 12569.5 - LD? (x, sin (sqrt(|x,))) 
Function 2.26 
Generalized 
Io Rastrigin's 
Function 


Sf) = - Wexp ( -0.2 (1/30L?_, x7)'7) - 
exp(1/30 X._, cos(2nx,)) + 20 +e 


: Griewank 3 
Fin | eee Sh) = EP, x2/4000 — []?,_, cos(x,/sqrt(i)) +1 ae | -600 <x, < 600 


Maximum number of function evaluation (FES) is set to 100.000 and the 
dimensions of the problems 30-D and 10-D for all the problems. 

Figure 2 depicts the convergence characteristics in terms of the best fitness 
value of the median run of each algorithm, which is plotted using a logarithmic 
scale to reduce the biggest and smallest values in the whole optimization process. 
Some of the illustrations have been omitted in order to save space. 


Sia Ackley’s Function 


6.2 Simulation Results and Performance Assessment 

To evaluate the performance of the algorithms each test functions were carried 
out 20 independent runs by the algorithms. After these runs the mean of the 
results and the standard deviation (STDEV) are recorded for each algorithms 
and test functions which are reported in Table 2. The best results are typed in 
bold. 

In Table 2, the proposed BAHC has demonstrated a better ability of global 
searching for all the benchmark functions except for function f, and f,. BAHC 
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achieves notable improvement on functions f,, f,, f,, f,. and f,. Although BAHC 
outperformed other algorithms for the function f,, it appears that all the 
algorithms trapped into similar local optima. It can be seen from Table 2 that 
using Hill Climbing with Bat Algorithm significantly increased the precision of 
the optimization results. 

For the function f,, BAHC performed very poorly comparing with NABA, 
CLSBA, and BA. f, is Quartic Function with noise which has ill-behaved 
surface. The locations of its minima and maxima are very close, so using Hill 
Climbing cost too much time to BAHC. 

BAHC showed better performance than other algorithms on f, f,, and f,, 
when the functions are 30 dimensions, but it didn’t show same performance 
when the functions are 10 dimensions. Because BAHC loses too much time to 
reach the local minima at the beginning of the algorithm. 

Furthermore, Fig.2 shows the convergence process of NABA, ABA, CLSBA, 
BA, and BAHC, conducted on some of the test functions from all the groups. 
The simulation results show the fast converging ability of BAHC in comparison 
with BA for most of the benchmark functions. Especially if one of the bat is on 
basin which contains global optima, the Hill Climbing part of BAHC increases 
the convergence speed rapidly as it can be seen for the functions, f,, f,,. and J, . 

A close inspection of Table 2 reveals that, out of 11 functions, BAHC 
outperformed all other contestant algorithms on 9 functions in terms of Mean. 
It can be clearly seen that BAHC can provide a better optimization solution for 
most of the functions encompassing different type of functions with high 
dimensions. Therefore, it can be said that BAHC has the best universality on 
different type of problems. Additionally, we can say that, NABA has the best 
performance in terms of the mean after BAHC. 


6.3 Algorithm Complexity 

This section describes the algorithm complexity of BA and BAHC as defined 
in (Liang et al., 2013). All the algorithms are developed in the Python 
environment and run ona PC with a 3.20 GHz CPU and 6.00 GB of RAM. For 
the random number generators, the Mersene Twister is used by Python in our 
simulations (Matsumoto and Nishimura 1998). Table 3 shows the BA and BAHC 
complexity on 10, 30, and 50 dimensions. As defined in (Liang et al., 2013), TO 
is the time calculated by running the following test problem: 


for i=1:1000000 
x=0.55+ (double) 1; 
X=X+X; 
x=x/2; 


X=x*xX; 
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X=sqrt(x); 

x=log(x); 

X=exp (Xx); 

x=x/(x+2); 
end 


T1 is the time to execute 200,000 evaluations of benchmark function f, by itself 
with D dimension, and T2 is the time to execute algorithms with 200,000 evaluations 
of f, in D dimensions. T2‘ is the mean T2 values of 5 runs. Based on Table 3, 
both T1 and T2‘ scaled linearly with the number of dimensions, as shown by the 
linear growth of (T2‘-T1)/TO for both algorithms. However, BAHC stayed behind 
the BA in all dimensions. Hill Climbing is slowing down BAHC. 


Table 2. Performance comparison of Means and Standard Deviations for 
NABA, ABA, CLSBA, BA, and BAHC 


eg PP 
Mean (StdDev) Mean (StdDev) Mean (StdDev) Mean (StdDev) Mean (StdDev) 


3 


30 
30 
30 

0 
30 
30 

0 
30 


5274,23(436.74) 5621.54(379.73) 5549.87(1275.91) 6384.07(376.02) 4158.53(185,24) 


223.61(116.29) 335.94(14.90) 186,40(21.63) 350.39( 189.86) 54.07(5.67) 
0.0045(0.0014) 5.8277(0.7727) 0.0047(0.0010) 0.0086(0.0028) 1.0357(0.5603) 


5.149(0.568) 8.546(0.853) 15.964(2.541) 12.347(0.982) 2.154(0.0254) 
48.69(5.8083) 58.12(6.2412) 52.82(8.1370) 54.09( 1.4475) 10.10(0.8899) 
0.0005(0.0002) 0.0817(0.0231) 0.0007(0,0005) 0.0009(0.0003) 0,001 1(0,0008) 


47.37(20.89) 54.79(11.76) 53.79(17.52) 39.79 7.233) 38.87(6.466) 
196.76(35.25) 317.21(35.27) 214.28(8.80) 246.92(1 1.27) 172.61(20.22) 
8.322(8.430) §.447(7.257) 15.52(6.904) 8.958(8.975) 15.34(6.827) 


1405.62(265.21) 1229.85(347.98) 1686.09(389.28) 1509.57(467.18) 1018.62(159.54) 


30 2 

30 19.81(0.2147) 19.96(8.1E-04) 20.00(0.2873) 19.94(0.0161) 19.72(0.1157) 
0.6892(0.0792) 0.8576(0.0565) 0.8168(0.0308) 0.7886(0.0581) 0.7518(0.0665) 

30 1.4284(0,0135) 1,4794(0.0979) 1,5004(0,0147) 1,4693(0.0067) 1.0791(0.0101) 


0.5905(0.0652) 0.3185(0.0431) 0.2924(0.0013) 2.3131(0.5116) 
63.36(68.99) 29.09(13.06) 12.26(1.111) 73.10(121.11) 8.844(0.874) 
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Table 3. Algorithms complexity 
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Fig. 2: Convergence characteristics of NABA, ABA, CLSBA, BA, and BAHC on the 
functions f,, ff, Sis Figs and f,,. 
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7. CONCLUSIONS AND FUTURE WORKS 


The present paper provides some information about BAHC, a hybrid of bat 
algorithm and hill climbing to solve high dimensional function optimization 
problems. Its performance was compared with those of NABA, ABA, CLSBA, 
and BA over 11 optimization problems. Results reported that the mixing BA 
with HC enhanced the optimization accuracy and improved the computation 
time of BA. Evaluation results obtained on optimization functions proved the 
effectiveness of the hybrid algorithm. 

Although the performance of BAHC outperforms BA and other algorithms 
in terms of accuracy and computation time for high-dimensional problems, 
BAHC does not perform so effectively for low dimensional problems which 
need low computation load. However BAHC spends too much time during hill 
climbing and then starts converging. As a future work convergence ability of 
BAHC for low-dimensional optimization problems would be studied. 
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ABSTRACT 


Rebuilding applications by reusing the existing ones is of primary importance for the 
IT developers. Every system uses resources such as Web pages, business information, 
etc. that can be represented in a computer-based system because the main purpose of 
IT department is making these resources accessible to its clients. Service architects 
and developers are interested in creating services that are implementable, maintainable, 
extensible and scalable. It is the RESTful design that promises this and a lot more. It 
has been more than a decade since the introduction of Representational State Transfer 
(REST), and it seems to become one of the most important technologies for Web 
applications. Every major development language now includes frameworks for building 
RESTful Web services. Considering the importance of REST while reusing services, 
we have conducted a survey in IT companies in Macedonia, which relies upon collection 
of empirical data. The survey aimed at exploring developers’ experience on REST 
Architecture in IT Companies and the data are here reported and evaluated. 
Keywords: REST Architecture, service reuse, resources, empirical data. 


1. INTRODUCTION 


The products or services offered by businesses today rely on information 
products developed by enterprise IT department. Time to market of a product 
is very significant for enterprises in order to be competitive in the market. Asa 
result, minimizing the production time by IT departments for a given product is 
essential for the enterprises (Algermissen, 2011). 
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The primary attraction for developers is to rebuild applications by reusing the 
existing ones. Reusing existing objects/components or making use of class 
libraries, functions, data, etc is as old as retelling a story (Imeri and Antovski, 
2012). To avoid duplicate designs of similar processes, abstract business 
processes were used to support for reusability. 

Every system uses resources, such as Web pages, business information, or 
anything that can be represented in a computer-based system. The purpose of 
a service is to provide to its clients an easy access to these resources. Also 
service architects and developers are interested in creating services that are 
easy to implement, maintain and are extensible, and scalable. It is the RESTful 
design that promises this and a lot more. 

It has been more than a decade since REST introduction, and it seem to 
become one of the most important technologies for Web applications. Its 
importance seems to grow very quickly since all technologies are moving 
towards Application Programming Interface (API) orientation. Every major 
development language now includes frameworks for building RESTful Web 
services. As such, it is important for Web developers and architects to have a 
clear understanding of Representational State Transfer (REST) and RESTful 
services. Many books and tools have been created but there is still a general 
lack of understanding its fundamentals as an architecture style. The reason 
perhaps could be found in the fact that REST was presented in a doctoral 
dissertation, with relatively high entry barriers for its understanding, or because 
the description used models that were more oriented towards documentation 
than to working practitioners (Wilde and Pautasso, 2011). 

Seeing the importance of REST while reusing services, we have conducted 
a survey in IT companies in Macedonia, which relies upon collection of empirical 
data. The survey aimed to explore developers’ experience on REST 
Architecture, its benefits and challenges from their viewpoint and see the 
possible contributors towards successful reuse, which could help increase of 
the knowledge and understanding of REST. We have prepared a questionnaire 
and sent it electronically to developers in 30 companies in Macedonia and 
other Balkan countries, but only the companies in Macedonia replayed to the 
questionnaire. Consequently, comparing the results obtained from IT companies 
in Macedonia and the region failed. 

The forthcoming part of the paper reports on the state of the art of the REST 
Architecture, the questionnaire and research questions and empirical results. 
In the end, the conclusions are drawn. 
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2. REST ARCHITECTURE 


Restful Services as architectural style lighter than SOAP-based Web Services, 
due to its simplicity, heterogeneity and web-based format, entities/resources 
are identified by unique URLs, including how resource states are addressed 
and transferred over HTTP by a wide range of clients written in different 
languages (Memeti et. al., 2015). 

REST originated at the intersection of academia and software development, 
among the architects of the World Wide Web (Adamczyk et al., 2011). 
Unfortunately, researchers have only recently started to work on RESTful 
services. No paper about RESTful Web services has been reported up to 2007. 
Information about RESTful services was reported only in 2010. 

The REST architectural style has been chosen for the resource centric 
approach because it could be easily applied and naturally mapped to it. The 
HTTP methods can be mapped to the Create (POST), Retrieve (GET), Update 
(PUT) and Delete (DELETE) operations and the XML documents can be 
used to provide a uniform resource representation (Szepielak, 2006). 

Many frameworks and tools for building RESTful Web services are available 
today. They are written in different programming languages and range from 
simple to quite sophisticated in their support of HTTP and other Web 
technologies. As they continue to improve, misunderstandings and violations 
present in today’s Web services will likely lessen (Wilde and Puatasso, 2011). 

When multiple providers offer the same service, a client has a choice and 
can select the most suitable one. Often, this choice comes down to the Quality 
of Service (QoS) parameters (Wilde, 2011). RESTful Web services today ignore 
QoS requirements; their only concern is providing functional interfaces. To 
add QoS parameters to RESTful services, a language for describing the 
parameters and a mechanism to incorporate the description in the HTTP payload 
is needed. Defining a standard QoS description language might benefit from 
the work in Semantic Web. Semantic Web ontologies define standard ways of 
interpreting information, such as QoS parameters, enabling all clients to interpret 
them the same way (Wilde, 2011). 

Building an integrated software environment in an enterprise often requires 
developing large amounts of Web services. The integration efforts can be greatly 
reduced by using a specialized framework for their development. Providing 
such tools that simplify software development in integration projects is essential 
for optimizing their efficiency and cost. One of the most important choices to 
make when building an integration solution is to select an appropriate integration 
approach and suitable technologies for its realization (Wilde, 2011). 

According to Fielding, (2000) the four principles of REST (called constraints) 
are: i) identification of resources, ii) manipulation of resources through 
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representations, tii) self-descriptive messages and, iv) hypermedia as the engine 
of application state. 

These principles describe the architecture of systems and interactions that 
make up the Web. The building blocks of the Web are called resources. A 
resource is anything that can be named as a target of hypertext (e.g., a file, a 
script, a collection of resources). In response to a request for a resource, the 
client receives a representation of that resource, which may have a different 
format than the resource owned by the server. Resources are manipulated via 
messages that have standard meanings; on the Web, these messages are the 
HTTP methods (Fielding, 2000). 

Activities that companies are supporting such as knowledge management 
and communication or decision making, the development approach such as 
object -oriented or component-based approach which is characterized by 
reusability (elements for re-used in other workflows), substitutability (alternative 
implementations easily inserted with precisely specified interfaces, and with 
the ability to verify and validate substitutions), extensibility and scalability (the 
ability to extend system component and to scale it, by increasing capabilities of 
individual components, with new functionalities), customizability (ability to add 
new features based on the needs of a particular domain), and composability 
(easy construction of complex solutions using basic components and the ability 
to estimate the efforts of system integration and compose) are of great 
importance for the developers. 

According to (Choi and Kim, 2008) the advantages of REST architecture 
are: i) availability: utilization of service which increases the reusability service, 
li) generality: a more generic service, since a reusable service could be 
implementable for unknown future requirements, iii) understandability: service 
and its interface are understandable. The more understandable it is, the more it 
is used, iv) functionality: to have a highest degree of guaranteed reuse, the 
service should pack with a complete range of functionality, v) reliability: since 
the service is to be used by many consumers, it should be reliable, be able to 
operate without any failure continuously and, vi) portability: the portable 
services which are to be consumed through any kind of environment without 
being modified increases the reusability service. 


3. EMPIRICAL ANALYSIS 


The present paper reports on usability of the REST architecture, its advantages 
and challenges. Consequently, a questionnaire with predefined answers and 
some open questions for general reflection was compiled. The research involved 
both qualitative and quantitative studies. In order to have clear understanding, 
at some points, some of the data were transformed from qualitative into 
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quantitative presentations. 

As regarding the sampling frame and size of our survey, the data presented 
in this paper are answers taken only from 10 software developing companies 
in Macedonia. The questionnaire was electronically sent to 20 companies in 
the region, but 10 companies in Macedonia replayed to the questions, only. 

Reported in Table 1, the questionnaire consists of 13 questions which have 
not been made before. Below are presented research questions of our survey. 

I. RQI. Do companies use REST Architecture? 

2. RQ2. Does developers’ experience with REST influence to the 
percentage of service reuse? 

3. RQ3. Which factors influence toward use of REST Architecture? 

4. RQ4. Which challenges and benefits are reported from companies? 


Table 1. The questionnaire 


Q1.Which Architectural Styles and Patterns are used by your organization? 

Q2.What is your experience with REST Architecture? 

Q3.Which activities are supporting the architectural styles used by your company? 

Q4.Which development approach is used by your organization? 

Q5.When a development team set up a project, would they start to go through old 
services or just start to write the new program application? 

Q6.What is the percentage of service reuse? 

Q7.How is the software production process carried out in your organization? 

Q8.While enhancing performances of existing applications, do you have any problems 
with permissions to use same data? 

Q9.If yes, how do you resolve the problem? 

Q10.Does the quality and flexibility of applications/services increased using REST 
(based on your experience)? 

Q11.Does the productivity of applications/services increases using REST (based on 
your experience)? 

Q12.Which advantages of REST Architecture are more important from your point of 
view? 

Q13.Which challenges do you face while using REST Architecture? 


3.1 Results 


The survey involved software development companies, with more than 3 
years working experience, whose application domain is information systems 
such as administrative, commerce and educational software and the data of 
the survey are in forthcoming part of the paper provided. 
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A) RQI. Do companies use REST Architecture? 

Answers to questions Q1 gave us information about the way companies 
develop software, i.e. which architectural style are using in their company. 
From the answers we could see that companies use different styles. Some 
have chosen one and some have chosen more than one style. The results were 
summed up and presented as percentage data to find the mostly used 
architecture. Depicted in the Figure 1, the results reported that 60% of them 
use REST architecture; 30% component-based applications; 20% data-centric; 
10% event/driven; 40% service-oriented and 60% client/server which gives 
the same percentage like REST usage. 


Service-oriented 
REST 


Event/driven 


Component-based 


Client/Server 


Percentage (%) 


Fig. 1: Percentage Architectural Styles used by Software Companies. 


B) RQ2. Does developers’ experience with REST influence to the 
percentage of service reuse? 


Since REST is relatively a new approach, our second research questions 
tries to see the relations among developers experience with REST and the 
percentage of service reuse. Based on surveys on software reuse (Imeri and 
Antovski, 2011) which is closely related to component based software reuse, 
we hypothetically concluded that the more experienced developers are, the 
higher the reuse is. The results obtained from the question Q2 reported that 
80% of the developers had high experience with REST. 20% of the developers 
had low experience. The results on reuse service obtained from question Q6 
reported different values (Figure 2). So, the developers’ experience with REST 
merely has a considerable impact on service reuse. 
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Fig. 2: REST experience vs. service reuse 


C) RQ3. Which factors influence toward use of REST Architecture? 


Information about the factors influencing the use of REST Architecture is 
from the questions Q3, Q4, Q5 and Q7 provided information. 

80% of companies are using Knowledge Management and Communication 
(KMC), and 20% of them Design Reasoning and Decision Making (DRDM). 

Regarding the development approach, 70% of companies have chosen object 
oriented (OO) paradigm and 30% of companies have chosen component-oriented 
(CO) paradigm. When a development team set up a new project, 60% of 
companies responded that while they develop new applications they go through 
old services (OS) whereas 40% start to write new program (NP). 

Finally, questions regarding the software production process, 50% of 
responders answered that the organization produce specific software (isolated 
products - IO) whereas 50% product families (PF), meaning that organization 
develops product that evolves over time, being somehow adopted for each 
new project (in-house products). 

Contingency tables (Table 2) were used to investigate the influence/dependence 
of the aforementioned factors by displaying the contingency distribution of 
variables. As we had 2x2 tables for the Chi-squared test, the Yates’ continuity 
correction was applied. 

We have applying the Yates’ Chi Squared Test (as one of statistical test 
applied to sets of categorical data), since the effect of Yates’ correction is to 
prevent overestimation of statistical significance for small data (https:// 
en. wikipedia.org/wiki/Yates%27s_correction_for_continuity, 2015), based on 
the fact that the number of our responders is small. 
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Table 2 reports the Contingency tables of observed and expected frequencies 
of different variables which influenced the use of REST Architecture such as: 
i) the interrelation between the uses of REST Architecture vs. Companies 
development approach, ii) the use of REST Architecture vs. Activities supporting 
the architectural style, iii) the use of REST Architecture vs. Software production 
process and, iv) the use of REST Architecture vs. Developing process. Lastly, 
the significance of the difference between the two variables is assessed with a 
variety of statistical tests including Yates’ chi-squared test. 

The critical value of S«? at a 5% significance level for 1 degree of freedom 
is seen to be 3,841 (https://en.wikipedia.org/wiki/Chi-squared_distribution). 
Below, are presented the calculated values of S«?. Almost in all our cases 
taken as factors that influence the use of REST Architecture, the calculated 
S«? value were smaller than the critical value of S«’, i.e., independent variables 
except the software production process found near the critical value of S«?, 
and can be taken as a dependent variable, such a factor that influence the use 
of REST Architecture. 


Table2. Contingency tables of factors influencing the use of REST. 


CO OO Total KMC DRDM_ Total 
Observed Others 4 Observed Others 4 
Total 3 sf 10 Total i 3 10 
CO OO Total KMC DRDM_ Total 
Expected REST ee 6 Expected REST al 6 
Others 4 Others 4 
Total 3 i 10 Total T i} 10 
w& Yetes’ 0,607143  Yetes’ 7 1,638889 
OS NP _ Total 10 PF Total 
REST 3 3 6 REST [4599255 « 
Observed Others eared 4 Observed Others ae i al 4 
Total 6 4 10 Total B) 5 10 
OS NP _ Total 10 PF Total 
Expected REST eye 6 Expected REST a ae | 6 
Others te | 4 Others | 4 
Total 6 4 10 Total 5 5 10 
eg Yetes’ 1,100694 ® Yetes’ a5 


D) RQ4. Which challenges and benefits are reported from companies? 


The challenges of REST consist of data permissions. Results reported that 
40% of developers have had problems with data permission due to their 
sensibility. 60 % of developers have not had problems with permission. Results 
are in Figure 3 depicted. 
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Fig. 3: Data permissions problems. 


Figure 4 depicts some of the characteristics of the REST architecture. 80% 
of the companies reported that they use the Rest architecture because of its 
flexibility. 60% of the companies reported that they use the Rest architecture 
because of its scalability. 
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Fig.4: REST Architecture benefits. 


4. CONCLUSIONS 


The presented investigation addresses the companies’ efficiency when meeting 
the customers’ requests and demands for new products. In general the empirical 
data collected represent an overview of the current situation among companies 
using the REST architecture. The most evident patterns that could be identified 
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from our analysis, such factors that influence toward the use of REST 
Architecture are presented in Table 2, where variable dependencies are 
calculated. 

The data statically (empirically) reported that the use of REST Architecture 
in Macedonian companies is relatively high, about 60%. Developers’ experience 
merely has any impact on service reuse. 

The factors that influence the use of REST architecture are Table 2 reported, 
showing that low probabilities of the tested variables dependence among factors 
we have chosen. A Yetes’ value (S«?= 1, 638889) is the calculated value among 
the use of REST architecture and activities supporting the architectural styles. 
As far as the development approach was concern, the calculated variables 
show that (S«*= 0, 607143) is the Yetes’ value among the development approach 
and the use of REST Architecture. In addition, a value (S«’= 1, 100694) is 
reported as a Yetes’ value while calculating the development approach such 
working through old services or write new program application as a activity 
supported by investigated companies, and lastly a Yetes’ value (S«= 3, 5) 
based on the software production process. Consequently, all chosen variables 
are independent from the use of REST architecture. The process of software 
production taken as a variable in our case influences the use of REST because 
of the low probability of dependence, considering that the calculated S«? value 
is near the critical value of S«’. 

The advantages of REST consist of flexibility and productivity. Results showed 
that about 80% of applications that involved REST architecture showed positive 
impact, which is in line with the principles of Fielding (2000). 
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ABSTRACT 


The present paper provides information about the properties of geotextiles based on 
their application area. Geotextiles are planar structures consisting of synthetic or natural 
polymeric material, used in contact with soil and other materials. Geotextiles are widely 
used in the construction industry due to their physical mechanical and hydraulic 
properties and durability. The present paper reports on the most important mechanical 
and hydraulic properties of geotextile used in specific applications areas along with 
their respective test methods based on Albanian and European Standards. 

Keywords: geotextile, properties, test methods, standards 


1. INTRODUCTION 


The properties characterizing the geotextiles make them fit for application in 
civil engineering such as filtration, drainage, reinforcement, separation, 
protection and fluid barrier or strength release purposes. Each application requires 
geotextiles performing one or more specific functions and meeting some specific 
criteria related to fibers’ characteristics and fabric construction, referred to as 
the structure of the geotextile (Karcz and Holtz 1988). 

In addition to endurance and durability, geotextile are mentioned for their 
physical, mechanical, hydraulic properties. For each property, geotextiles are 
submitted to the respective testing method. Geotextiles is the combination of 
the physical form of the polymer fibers, their textile construction and the polymer 
chemical characteristics. In addition to the aforementioned properties, geotextile 
testing might be characterized as index or performance as a general value 
from which a specific property could be qualitatively assessed and used for 
general characterization of a geotextile product could be provided. Moreover, 
these tests do not provide data that can be directly used for design purposes 
and are performed on the geotextile alone, or in-isolation. For instant assessment 


126 Ilda KOLA, Blerina KOLGJINI 2015/Vol. XX (2) 


geotextiles are required to be investigated along with associating material such 
as soil, etc. As a result, information about the expected behavior of a geotextile 
in an engineered system could be obtained (Holtz et al., 1998). 

Physical properties help characterize the geotextile in the as-received, 
manufactured condition. Data on the physical index could be obtained via index 
testing. Common physical properties are specific gravity, mass per unit area, 
thickness and stiffness. 

Long-term performance or durability of geotextiles is also affected by 
geotextile degradation via different mechanisms including ultraviolet light 
(sunlight), chemical reactions with geotextile polymers, and/or thermal 
degradation. Degradation testing is important in determining the ultimate lifetime 
of a geotextile in an engineered system (Zornberg and Thompson, 2012). 

Mechanical responses and filtration ability are the most important properties 
of geotextile as they are a means to address working efficiency (Karcz and 
Holtz 1988). The present paper provides information about these properties. 


2. Mechanical properties 


The mechanical response of a geotextile depends on the orientation and 
regularity of the fibres and the type of polymer from which it is made. As 
aforementioned said, the geotextile must perform work in a stressed environment 
and be able to resist damage in an arduous environment. 

Usually the stressed environment is known and the textile is selected on the 
basis of some criteria to cope with the expected imposed stresses and its ability 
to absorb those stresses over the proposed lifetime of the structure without 
straining more than a predetermined amount. 

Figure 1 compares the tensile behaviors of a range of geotextiles (Horrocks 
and Anand, 2000). On the other hand, damage can be caused on site when 
constructing (e.g. accidental tracking from vehicles) or in-situ, when using 
geotextile (e.g. punching through geotextiles by overlying angular stone). 

Performance ability is based on the textile stiffness in tension and its ability 
to resist creep failure under any given load condition. The ability to resist damage 
is complex, being a function of the fibre’s ability to resist rupture and the 
construction of the fabric, which determines how stresses may be concentrated 
and relieved (Horrocks and Anand, 2000). The weight or area density of the 
fabric is an indicator of mechanical performance only within specific groups of 
textiles, but not between one type of construction and another. For example, 
within the range of needle punched continuous filament polyester fabrics, weight 
will correlate with tensile stiffness, while a woven fabric with a given area 
density will be much stiffer than an equivalent weight needle punched structure. 
So, the construction controls the performance. Therefore, weight cannot be 
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used alone as a criterion in specifying textiles for civil engineering use, but in 
combination with other specified factors, weight can be a useful indication of 
the kind of product required for a particular purpose (Christopher and Holtz, 
1985). 


Force (kN m™ width) 


o © 2 OD 0 6 
Strain % 


Fig.1: Typical ultimate stress-strain failure levels (a) of high strength and (b) of medium 
strength PES woven geotextiles used for embankment support and soil reinforcement, (c) of 
geogrids and lower strength polyester woven geotextiles used for soil reinforcement and (d) 

of low strength, highly extensible nonwoven geotextiles used for separation and filtration 


(Horrocks and Anand, 2000). 


One of the most important strength properties for a geotextile is the tensile 
strength. The breaking strength of a standard width of fabric or ‘ultimate strip 
tensile failure strength’ is universally quoted in the manufacturers’ literature to 
describe the ‘strength’ of their textiles. But this is of very limited use in terms 
of design. Rather, a strength at a given small strain level will be the design 
requirement. Therefore, the tensile resistance or modulus of the textile at say, 
2%, 4%, and 6% strain is much more valuable. To understand the load-strain 
characteristic, it is important to consider the complete stress-strain curve. It is 
also important to consider the nature of the test and the testing environment 
(Holtz et al., 1998). Ideally, continuous stress—strain curves should be provided 
for engineers, to enable them to design stress resisting structures properly. 
Geotextile tensile strength is measured by clamping two opposite ends of a 
geotextile specimen in a mechanical testing machine and stretching the specimen 
until failure occurs. Typically, both the force applied to the geotextile and the 
geotextile strain are measured allowing observation of the stress-strain curve 
and development of associated modulus (Zornberg and Thompson, 2012). 

Stress—strain curves, as shown in Fig. 1 and 2, may well comprise a high 
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strain sector, contributed by the textile structure straightening out, and a low 
strain sector, contributed by the straightened polymer taking the stress. In Figure 
2, geogrids (a) absorb the imposed stresses immediately, giving a high initial 
modulus. Later, the curve flattens. Woven fabrics (b) exhibit initial straightening 
of warp fibres which produces a low initial modulus. Later the modulus increases 
as the straightened polymer fibres take the stress directly. Nonwovens (c) give 
a curvilinear curve, because extension is primarily resisted by straightening 
and realignment of the random fibre directions (Horrocks and Anand, 2000). 
Geotextiles might have different strengths in different directions. Therefore, 
tests should be conducted in both principal directions. The wide-width tensile 
test S SH EN ISO 10319 is the most commonly specified in Albanian and 
European Standards for testing geotextiles used in almost all the applications 
performing all functions (Drejtoria e Pérgjithshme e Standardizimit, 2008). 


a 


Force (kN m™ width) 


Strain % 


Fig.2: Different stress—strain curve shapes exhibited by the three main types of 
geosynthetic construction. (a) Geogrids, (b) Woven fabrics, (c) Nonwovens. 


Tensile creep is a time-dependent mechanical property and strain at constant 
load. Creep can cause the physical failure of a geotextile if it is held under too 
high mechanical stress. Creep tests can be run for any of the tensile test types, 
but are most frequently performed on a wide strip specimen by applying a 
constant load for a sustained period. Short-term creep strain is strongly influenced 
by the geotextile structure. Woven geotextiles have the least; heat-bonded 
geotextiles have intermediate; and needled geotextiles have the most. Longer- 
term creep rates are controlled by structure and polymer type. The creep limit 
is the most important creep characteristic. It is the load per unit width above 
which the geotextile will creep to rupture. The creep limit is controlled by the 
polymer and for polyester it is approximately 60% of its ultimate strength, for 
polyethylene about 40% of its ultimate strength and for polypropylene around 
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20% of its ultimate strength. Therefore, for example, a polyester fabric with an 
ultimate tensile strength of 1OOkNm-1 width cannot be loaded under a long 
term stress of more than 60 kNm-1. The higher the level of imposed stress 
above this point, the more rapid will be the onset of creep failure. Figure 3 
depicts the safe loading limits for most commonly used geotextiles (Horrocks 
and Anand, 2000). The S SH ISO 13431 ‘Determination of tensile creep 
rupture behavior’ is defined in the Albanian Standards as a means to address 
tensile creep for a number of specific applications (Drejtoria e Pérgjithshme e 
Standardizimit, 2008). 
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Fig. 3: Approximate limits of creep resistance for different geosynthetic polymer 
constructions (Horrocks and Anand. 2000). 


Other tests for geotextile strength are focused on ensuring that a geotextile 
is strong enough to survive installation stresses, often the most severe stresses 
placed on the geotextile over its lifetime. The most commonly specified tests 
cover tear strength and puncture strength (Zornberg and Thompson, 2012). 

Often, due to manufacturing size limitations, it is desirable to join multiple 
geotextile sections together when placed in an engineered system. The most 
common method of joining geotextiles is by sewing them. If sewing is used 
when installing geotextiles, characterization of the strength of a given seam 
type would be important for ensuring the integrity of the sewn sections. This is 
carried out by testing a tensile strength on a sewn geotextile sample as defined 
by the Albanian Standard, S SH EN ISO 10321 (Drejtoria e Pérgjithshme e 
Standardizimit, 2008). 
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3. Hydraulic Properties 


A geotextile is similar to a soil in that it has voids (pores) and particles (filaments 
and fibers). However, the geometric relationship between filaments and voids 
is more complex than in soils due to the shape and arrangement of the filaments 
and the compressibility of the structure with geotextiles. Since pore size can be 
directly measured, relatively simple relationships between the pore sizes and 
particle sizes of the soil to be retained can be developed. Three simple filtration 
concepts are used in the design process: 1) if the size of the largest pore in the 
geotextile filter is smaller than the larger particles of soil, the soil will be retained 
by the filter. As with graded granular filters, the larger particles of soil will form 
a filter bridge over the hole, which in turn, filters smaller particles of soil, which 
then retain the soil and prevent piping (Figure 4) (Holtz et al., 1998), ii) if the 
smaller openings in the geotextile are sufficiently large enough to allow smaller 
particles of soil to pass through the filter, then the geotextile will not blind or 
clog (Figure 5) (Holtz et al., 1998) and, iii) a large number of openings should 
be present in the geotextile, so proper flow can be maintained even if some of 
the openings later become plugged (Holtz et al., 1998). 

The following design criteria for geotextiles must be met: i) soil retention, ii) 
water permeability and, iii) clogging resistance. 

To perform effectively, the geotextile must also survive the installation process 
(survivability criterion) (Zhou, 1998). 
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Fig. 4 Filter bridge formation Fig. 5 Clogging and blinding. 


Hydraulic properties are related to the pore size distribution of the geotextile 
and correspondingly its ability to retain soil particles over the life of the project 
while allowing water to pass (Lee and Bourdeau, 2006). 

The ability of water to pass through a geotextile is determined from its hydraulic 
conductivity (coefficient of permeability, k), as measured in a permeability 
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test. The permeability of geotextiles can vary immensely, depending upon the 
construction of the fabric. National and international standards have been set 
up for the measurement of permeability that is required, most often at right 
angles to the plane of the textile (crossflow), but also along the plane of the 
textile (in-plane flow, called transmissivity). Albanian Standards specifying 
those tests are S SH EN ISO 11058 — Water permeability normal to the plane 
and S SH EN ISO 12958 — Waterflow capacity in the plane (Drejtoria e 
Pérgjithshme e Standardizimit, 2008). 

Water flowing is very important for civil engineering earthworks as it prevents 
the pressure buildup of water. Coefficient permeability is defined as the product 
of permeability considering its dimension in the flow direction. The units are 
nationalized in metres per second. Effectively the coefficient is a velocity, 
indicating the flow velocity of the water through the textile. 

Due to the compressibility of geotextiles, engineers also use a coefficient 
called the permittivity, which defines the theoretical permeability irrespective 
of the thickness of the fabric. The permittivity, w (permeability divided by 
thickness), is often determined from the test and used to directly evaluate flow 
capacity. Permittivity is defined as the cross-plane permeability of a geotextile. 
In other words, permittivity is a measure of the ease at which water may flow 
through the geotextile. Permittivity is obtained with S SH EN ISO 11058, by 
measuring the flow of water, under a given head, moving perpendicularly through 
a geotextile (Drejtoria e Pérgjithshme e Standardizimit. 2008). Values are 
reported in units of sec’', allowing for a more traditional value of permeability 
to be obtained by multiplying the permittivity by the thickness of the geotextile. 

Permeability criteria are established based on the critical nature and severity 
of the application. For noncritical, less severe applications, the permeability of 
the fabric (k,) is required to be at least greater than the permeability of the soil 
(k,) the fabric is to retain, k, >k,. For critical-severe applications the fabric 
permeability must be at least ten times greater than the permeability of the soil, 
k, > 10k, (Karez and Holtz 1988). 

Permeability criteria for woven geotextiles are in terms of the percent open 
area (POA). When the protected soil has less than 5 percent passing the No. 
200 sieve, the POA should be equal to or greater than 10 percent. When the 
protected soil has more than 5 percent but less than 85 percent passing the No. 
200 sieve, the POA should be equal to or greater than 4 percent (Joint 
Departments of the Army and Air Force, 1995). 

The ability of a geotextile to retain soil particles is directly related to its 
apparent opening size (AOS) which is the apparent largest hole in the 
geotextile. The AOS value is equal to the size of the largest particle that can 
effectively pass through the geotextile in a dry sieving test. It is also called 
filtration opening size (FOS). 
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The procedure for matching a textile to the soil, in order to achieve stability 
under difficult hydraulic conditions, is to use a textile whose largest holes are 
equal in diameter to the largest particles of the soil (O,, = D,,). Where hydraulic 
conditions are less demanding, the diameter of the largest textile holes can be 
up to five times larger than the largest soil particles (O,, = 5D,,). The relationship 
between O,, and D,, is in figure 6 depicted. Particularly difficult hydraulic 
conditions exist in the soil: 1) when under wave attack, ii) where the soil is 
loosely packed (low bulk density), iii) where the soil is of uniform particle size, 
or iv) where the hydraulic gradients are high. Lack of these features defines 
undemanding conditions. Between the two extremes lies a continuum of variation 
which requires experience and judgment in the specification of the appropriate 
O,, size for any given application (Horrocks and Anand 2000). 


Geotextile 
90% of all particles in the 
soil have a diameter 
smaller than this particle, 
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90% of all holes in the Jt 
textile are smaller than | ) 
this hole, so it is called | i) 


the Ovo 


Fig. 6: Relationship between O,, and D,, 


The largest hole sizes and largest particle sizes are assessed by consideration 
of the largest elements of the fabric and soil. Measuring the largest particles of 
a soil is achieved by passing the soil through standard sieves. In order to assess 
a realistic indication of the larger particle diameters, a notional size is adopted 
of the sieve size through which 90% of the soil passes. This dimension is known 
as the D,, by convention. Similarly, an indication of the largest holes in a textile 
is taken as the 90% of the biggest holes in the fabric, the O,,. Even under ideal 
conditions, if the O,, pore size is bigger than 5D,,, then so called piping will 
take place. The textile O,, pore size should be reduced from 5D,, towards D,, 
as the ground and hydraulic conditions deteriorate. The current standard in 
Albania for measuring the characteristic opening size is SSH EN ISO 12956 
(Drejtoria e Pérgjithshme e Standardizimit. 2008). 
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4. Situation in Albania 


The interest on geotextiles has been lately worldwide increased, especially in 
developing countries, including Albania. Some applications of geotexiles in Albania 
are: in the Levan -Dames road (Maccaferri Corporate, 2012), where a non- 
woven geotextile was placed on top of the pile caps to protect the gravel drain 
from becoming contaminated by fine soils transported within the ground water. 
Also in the Gjirokastér -Tepelené road (Maccaferri Corporate, 2003), where the 
geotextile is placed at the bottom of the structure in order to save it from the river 
flow erosion. Moreover, along the cut of the slope was installed a drainage 
composite to guarantee the drainage of the rainfalls. Both segments are part of 
the North South Corridor in Albania between Montenegro in the North and Greece 
in the south. Another application is in the Durrés - Kukés - Prishtiné Motorway, 
part of the South — East Europe Core Road Network, where the geotextile is 
placed at the bottom of the embankment to stabilize soils (Gjetvaj et al. 2012). 

Short and long terms properties of geotextiles are very important for civil 
engineering and the use of geotextile should be based on standard for every 
application. Testing geotextiles and evaluating their applicability would be very 
important for the entire structure where they are applied. 

Currently, European standards have been adapted by the Albanian institutions 
involved in the area. However, their implementation has still to be defined as 
reviewing the Technical Construction Legislation with regard to some specifics 
about geotextiles and their use is necessary. Once reviewing has been made, 
including geotextiles in the list of construction materials would be appropriate. 
Here, evaluation of conformity is required in a normative way. 

Table 1 reports on the standard test methods used to consider the properties here 
reported as defined by the Albanian Standards based on the European standards. 


Table 1. Geotextile properties and their associated S SH EN ISO Standard 
Test Methods. 


Characteristic Test method Functions 


s Pr B and STR" 
Tensile strength SSH EN ISO 10319 x x x x x 
S SH EN ISO 10319 ) ) x x 
Elongation at maximum load 
S SH EN ISO 12236 > x x x 
Static pancture resistance (CBR test) 
S SH EN ISO 13433 c x x x 


Dynamic perforation resistance 


S SH EN ISO 11058 


Water permeability normal to the plane 


S SH EN ISO 12956 


Characteristic opening size 


S SH EN ISO 12958 x 


Water flow capacity in the plane 


* F — filtration; D — drainage; R — reinforcement; S — separation; P — protection, B and STR — fluid barrier and strength release 
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5. CONCLUSIONS 


Geotextiles properties make them fit for civil engineering use. Mechanical 
and hydraulic properties are the most important ones. 

Mechanical properties refer to the ability of a textile to perform work in a 
stressed environment and to resist creep failure under any given load condition. 
They provide appreciative geotextile strength and/or compressibility under 
varying loads. The most important are tensile strength and tensile creep. Load- 
strain and the complete stress-strain curve are closely related. Creep can cause 
the physical failure of a geotextile if it is held under too high a mechanical 
stress. 

Hydraulic properties relate to the ability of geotextile to retain soil particles 
over the life of the project while allowing water to pass. The ability of water to 
pass through a geotextile is related to characteristics such as permeability, 
permittivity, transmissivity, while retaining soil particles is related to the Apparent 
Opening Size of geotextiles and its relation to the particle soil size. 

The data here reported could be used to develop a properties classifying 
system in order to define geotextiles for routine applications in Albania. The 
classification system was carried out for the first time by the French Committee 
of Geotextiles and for the design of a number of routine geotextile applications. 
This system has been used by other countries to develop their own systems. 
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HEALTH THROUGH NUTRITION: FACTORS OF 
NUTRITION AND THEIR IMPACT ON MORBIDITY 


November 21, 2015 


Making Science work for the benefit of the Society 


The Section of Natural and Technical Sciences of the Albanian Academy of 
Sciences (hereafter reffered to as Academy) organized on November 21* in 
close partnership with the University of Medicine of Tirana, Ministry of Health, 
Institute of Public Health and Compulsory Health Insurance Fund, the Scientific 
Conference “Health through nutrition: factors of nutrition and their impact on 
morbidity”. As the Academy through its Section aims at introducing new 
initiatives in the area of research and education to suppor t current needs of 
the country and helping find solution to many crucial issues concerning the 
development of the country by providing high state institutions with appropriate 
advice and expertise, the conference run in line with the Government Program 
and aimed at: 1) evidencing the nutrition’s impact on morbidity of the population 
in Albania and prevention from specific diseases, ii) involving all the stakeholders 
and raising awareness on the role the food palys in morbidity reduction, iti) 
exchaning inforamtion, buiding and raising the capacities among the stakeholders 
involved in the area. 

The beneficary target groups were the Ministry of Health of Albania, 
Compulsory Health Insurance Fund, Regional Health Directories, Primary 
Service of Family Doctors, wide population and food business. 

At the end of the Conference, the Conference Declaration Chart was 
approved and recommandations to state institutions involved in the area were 
made. 
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ALBERT EINSTEIN AND THE GENERAL THEORY OF 
RELATIVITY-IN THE FRAMEWORK OF SCIENCE 
EDUCATION PROGRAM 


November 27, 2015 


2015 is called the year of the festival of light because it marks the 100" 
anniversary of the General Theory of Relativity by Albert Einstein, one of the 
most eminent personalities of all times. Consequently, the Section of Natural 
and Technical Sciences organized in the framework Science Education Program, 
the International Conference “Albert Einstein and the General Theory of 
Relativity”. 

He is the fuzzy haired genius who unlocked the secrets of the universe. He 
was also a great peacemaker. Thanks to Albert’s great E = mc? equation, the 
world has nuclear power, radiation therapy, and gamma ray scans. His 
explanation of the Photoelectric Effect has led to fluorescent lights, automatic 
doors, lasers, TV, computers and much more. Children, also have to be grateful 
to him any time they play video games. 

A list of guests included Mr. Dhori Kule, Rector of the University of Tirana, 
Albania, Ligor Nikolla, Dean of the Faculty of Ingeneerin mathmatics and 
Physics, Mynyr Koni, Dean of the Faculty of Natural and Technical Sciences, 
Academitians, professors of physics and students, members of the Albanian 
parliament and many representatives of higher education and research institutions 
and His Excellency Boaz Rodkin, Ambassador Extraordinary and Plenipotentiary 
of Israel in Albania. 

Of great interest were the papers of Prof. Halil Sykja, Prof. Dritan Spahiu, 
Prof. Dr. Mimoza Hafizi, Dr. Lorenzo Amati from the Italian Institute for 
Astrophysics who shed light on some of the achievements in the area of physics 
and astrophysics due to Albert Enstein’s theory. 

In the end, Dr Peter McGrath - Coordinator of the Science Diplomacy 
Program run by the World Academy of Sciences (TWAS) for the advancement 
of science in developing countries provided information about the legacy of 
Einstein as diplomat scientist. 
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IN VITRO CULTURE —- FROM SCIENTIFIC RESEARCH TO 
MARKET; THE NEED FOR A VIRTUAL NATIONAL 
NETWORK OF PLANT BIOTECHNOLOGY 


December 3, 2015 
Hall “‘Aleks Buda’, Albanian Academy of Sciences 


The Academy of Sciences is dedicated to promoting the exellence in science, 
its benefits to the society and providing scientific advice for policy. Consequently, 
the Section of Natural and Technical Sciences of the Albanian Academy of 
Sciences organised in close partnership with the Ministry of Agriculture, Rural 
Development and Water Administration; the Center of Agricultural Technology 
Transfer, Vlora; University of Tirana and Agricultural University of Tirana; 
Academy of Sciences and Arts of Kosovo and the University, Ukshin Hoti, 
Prizren, Kosovo the Scientific Symposium “Jn vitro culture — from scientific 
research to market; the needs for a Virtual National Network of Plant 
Biotechnology”. The Symposium aimed at araising awareness on the relationship 
between technologic development and products and services of high quality in 
line with Government Program. 

The main topics adressed were as following: 1) phytosanitary status of fruit 
trees in Albania, production of virus-free saplings applying the meristem culture 
technique and making certified material available to farmers, ii) current diagonstic 
opportunities for phytopatologies of cultivated plants applaying molecular 
methods, iii) application of tissue culture as an alternative method of medium 
and long term conservation of plant germplasm and a variability source for 
genetic improvement and breeding of cultivated plants, iv) use of in vitro 
technique as a means to address genetic plant improvement to create new 
plant forms of high quality and, v) improving and fostering biotechnological 
research in the area of in vitro culture in Albania by creating the Virtual National 
Network of Plant Biotechnology; financial support intented to be over-arching, 
in order to support the specific scientific areas would be of crucial importance. 

In the area of biotechnology, Albania is characterized by an excellent 
infrastructure and well-equipped research facilities with highly qualified staff: 
1) scientific laboratories of cell and tissue culture and molecular biotechnology, 
Faculty of Natural Sciences, University of Tirana, set up by budgetary funds. 
Up to 2007, these laboratories were assets of the Institute of Biological 
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Research, Albanian Academy of Sciences. Since then, these laboratories have 
been infrastructure for Bachelor, Masters and PhD programs and thesis in the 
area of in vitro technology of germplasm of spontaneous species which are an 
interesting component of Albanian flora (endemic, rare, threatened and with 
high economic interest species as the medicinal-aromatic plants, authocthonous 
fruit trees species etc.), identification of plants viruses, evalutation of genetic 
polymorphism in subspecies and interspecific level, and the evidence of 
somaklonal variations of in vitro plants. The skillfull staff make this Department 
very important for the University of Tirana, ii) the complex of scientific and 
commercial laboratories with greenhouses for production of fruit seedlings 
(former Institute of Pomology, Vlora, currently the Center of Agricultural 
Technology Transfer), which primary purpose is the certification of fruit species 
rootstocks (in the framework of Italy-Albania Project in the collaboration with 
CIHEAM (Mediterranean Agronomic Institute), Bari, Italy, in the framework 
program INTERREG III, 2001-2011, of the European Union (EU). These 
laboratories were set up with substantial funds to shore up the support for the 
faculty. Up to 20 thousand certified seedlings of fruit species rootstocks were 
produced up to 2011, when the production stopped, iii) laboratory of plant genetic 
bank, State Agency for Seeds and Saplings (Project of the World Bank in the 
collaboration with the Project of Agricultural Services, ASP) (2003 — 2007), 
which is part of the Plant Genetic Resources Center since 2007. Unfortuantely, 
the laboratory has not been put into function, yet, regardless considerable 
investments and the efforts made to put the laboratory into function., iv) scientific 
and didactic laboratory, Department of Agricultural Production, Faculty of 
Agriculture and Environment, Agriculture University of Tirana (TEMPUS 
Program) (2011-). Eventhough the laboratory is well-equipped, further training 
would be necessary for a better scientific research in the field of in vitro 
cultures and, v) didactic laboratory, Department of Horticulture and Landscape 
Architecture, Faculty of Agriculture and Environment, Agriculture University 
of Tirana (project of Ministry of Education and Sports, UNILAB), which 
additional equipment and specialized staff are of fundamental importance. 

In his welcome speech, the Head of the Section of Natural and Technical 
Sciences, Acad. Prof. Dr Salvatore Bushati emphasized the impact of 
technologic development on biotechnology and Science education Program. In 
addition to prof. Bushati, Prof. Dr. Ramé Vataj from the University, Ukshin 
Hoti, Prizren, Kosovo; Prof. Dr. Mynyr Koni Dean of the Faculty of Natural 
Sciences, University of Tirana; Prof. Dr. Ardian Maci by the Faculty of 
Agriculture and Environment, Agricultural University of Tirana and Mr. Flamur 
Hysi on behalf of the Ministry of Agriculture, Rural Development and Water 
Administration welcomed the participants to the conference. In this Session, 
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Acad. Ass. Efigjeni Kongjika held main presentation “State and perspective 
of scientific research in the field of plant biotechnology (in vitro culture) 
in the World and Albania”. 

After presenting the orchard phytosanitary situation and the risk of viral 
diseases spread in nursery in Europe and Albania by various specialists from 
several institutions (CERTIS Company, Italy; CIHEAM, Mediterranean 
Agronomic Institute, Bari, Italy; Faculty of Agriculture and Environment, Tirana 
Agricultural University and Faculty of Natural Sciences, Tirana University), 
the specialists of the Center of Agricultural Technology Transfer, Vlora etc. 
presented scientific data on the in vitro culture application for the production 
and micropropagation of sanitated plant material. Aspecial attention was paid 
to the presentation of these specialists with regard to the possibilities of 
introducing sanitated biotechnologic product in the Albanian market and to 
constructe new private commercial laboratories. 

Many specialists of the Academy of Sciences, University of Tirana, 
Agricultural University of Tirana and the University, Ukshin Hoti, Prizren 
reported their long experience in the field of in vitro cultures on the 
micropropagation, meristem culture and in vitro conservation of plant 
germplasm. The possibilities of using new in vitro techniques and the perspective 
of use of modern biotechnologies at the laboratories of molecular biotechnology 
and in vitro cultures of the University of Tirana and Agricultural University of 
Tirana. 

A group of specialists from different institutions emphasized the need for a 
National Virtual Network of Plant Biotechnology as a means to address 
contemporary research in the field of plant biotechnology. The network would 
play the role of an ad-hoc committee where some of the most eminent 
personalities would come together to discuss the long-term needs, and tendencies 
of biotechnology. 

At the end of the symposium, a chart with recommendations from the auditor 
to the scientific community and policymakers was approved. 
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ANDREA PIERONI AWARDED HONORARY MEMBER OF 
THEALBANIANACADEMY OF SCIENCES 


December 7, 2015 


The Section of Natural and Technical Sciences of the Albanian Academy 
of Sciences is honored and priviledged to announce Prof. Dr Andrea Pieroni 
Honorary Member of the Academy of Sciences of Albania. 

Prof. Andrea Pieroni, after receiving his masters in Pharmacy and trained 
in Medical Botany at the University of Pisa, earned his doctorate in Life Science, 
University of Bonn, Germany. A Research Assistant at the University of London 
since 2000, as well as Lecturer and Senior Lecturer at the School of Life 
Sciences of the University of Bradford in Northern England since 2003, he 
returned to Italy in January 2009. 

Currently, he is a Professor for the Environmental and Applied Botany 
subject area (food botany, ethnobotany, and ethnobiology), University of 
Gastronomic Sciences, Pollenzo/Bra, Italy. Professor Pieroni has been the 
President of the International Society of Ethnobiology and currently he is Editor- 
in-Chief of the Journal of Ethnobiology and Ethnomedicine. He is also a 
board member of several International Ethnoscientific Societies and publication 
editors, as: Journal of Ethnobiology and Ethnomedicine (editor-in-chief), 
Journal of Ethnopharmacology, Journal of Ethnobiology, Ethnobotany 
Research and Applications, Food and Foodways, Ethnobiology and 
Conservation. 

He has been the PI. of the first collaborative research project funded by 
the EU Commission focused on traditional plant knowledge in the Mediterranean 
(RUBIA, 2003-2006). His research focus is on gastronomic and medical 
ethnobotany in the Mediterranean area, the Balkans, and among migrant 
communities in Western Europe. 

His courses are focused on: i) Food Biodiversity Sciences (in English), 11) 
Ethnobotany and Ethnobiology (in English), iii) Plant Sciences (in the Italian 
language), iv) Molecular basis and physiology of taste and Nutraceuticals (in 
the Italian language) and, v) Ecology(in the Italian language). 

His main research topics are: 1) ethnobiology: cross-cultural food and medical 
ethnobotany, i1) Ethnomedicine: folk medical practices, iii) human ecology: local 
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ecological knowledge and its mechanisms of adaptation and, iv) ethnoveterinary 
among pastoralists. 

Studies areas/communities concern linguistic and/or religious minority groups 
and cultural boundaries in South and South-Eastern Europe and migrant 
communities in Western Europe 

Research areas concern ethnobotany, ethnobiology, ethnomedicine, 
ethnoveterinary, social anthropology, migration studies, ethno-ecology, human 
ecology, gastronomy, ethnozoology, Balkan studies, Albanian studies, folk 
medicine, local gastronomy, medicinal anthropology, food anthropology, linguistic 
minority, Rumania (anthropology), gastronomical literature, gender, culture, 
emmigration and diaspora, ethnobiological studies, ethnobiological, 
ethnobiological for Albania. 

Community Studies involve linguistics and religious minority groups, and 
their cultural border in South/South-East Europe and immigrant communities in 
Western Europe, as well. 

His research focus is on Ethnobotanic Medicine and in the Mediterranean, 
Balkan and Immigrant communities in Western Europe. Current studies center 
on plant use as medicine and food, from the linguistic minorities in Southern 
Italy, Balkan and in immigrant communities (Eastern Europenians, South Asian, 
Turkish etc.) in Europe. He was a leader and pioneer in the first research 
project financed from the European Commission, focused on traditional 
medicinal knowledge in the Mediterranean (RUBIA- Ethnobotanical and 
ethnographic inheritance of the Mediterranean traditional technology, tools and 
uses of neglected and wild plants for food, medicine, textile and coloring) 2003- 
2006. 
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PROF. DR NIKO PANO AWARDED 
“HONOR OF THE ACADEMY” 


December 22, 2015 


The Albanian Academy of Sciences based on the proposal from the Section 
of Natural and Technical Sciences awarded Prof. Dr Niko Pano “Honor of the 
Academy”. 

Niko Pano is an eminent personality in the area of hydrology, internationally 
recognized as well. His early work in hydrology and the sub-fields is related to 
water and energy economics. In addition, he is well known for his efforts to 
foster intellectual pursuits of the staff of the first Institute of Hydrometeorology 
in Albania in the area of scientific research. 

He is the author of the “Water resources of Albania’, published by the 
Albanian Academy of Sciences. 


